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CHAPTER  I 
INTRODUCTION 

Phase  Change  and  the  Critical  Point 

Phase  changes  are  common  occurrences  in  matter.    Examples  are  the 
melting  of  ice,  the  evaporation  of  water,  and  the  sublimation  of  COg. 
These  changes  occur  slowly  and  in  small  portions  of  the  substance.  On 
the  other  hand,  the  domain  magnetization  of  iron  slowly  decreases  as  the 
temperature  is  increased,  until,  at  a  temperature  of  1043  F>  tne  magne- 
tization suddenly  vanishes.    This  transition  is  due  to  the  disappearance 
of  long  range  ordering  of  the  spin  magnetic  moment  of  the  3d  electrons. 

The  transition  in  iron  is  an  example  of  a  phase  transition  brought 
about  by  the  mutual  cooperative  interaction  of  many  particles.  Other 
examples  include  the  critical  point  of  liquid-vapor  transitions,  the 
Neel  point  of  antiferromagnets,  and  the  lambda-point  in  liquid  helium. 
All  of  these  are  characterized  by  an  order  parameter  describing  the  co- 
operative effects  which  increases  rapidly  below  a  critical  temperature. 
Several  such  systems  are  listed  in  Table  I.  taken  from  Kadarioff,  et  al.J 
and  Huang. 

A  ferromagnet  is  characterized  by  the  parallel  alignment  of  the 
spins  of  the  magnetic  ions  below  a  critical  or  Curie  temperature,  T  . 
Thermal  agitation  randomizes  the  spins  at  temperatures  above  the  Curie 
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Table  1.    Cooperative  Systems  and  the  Corresponding 

Order  Parameters 


System 

liquid  -  gas3 
ferromagnet 
antiferromagnet 
superconduc  tor 
super f luida 
ferroelectric3 

a 

phase  separation 
binary  alloy  (bcc) 


Order  Parameter 
density 

domain  magnetization 

sublattice  magnetization 

complex  gap  parameter 

condensate  wave  function 

lattice  polarization 

concentration 

2x  (number  of  similar 

nearest  neighbors)  + 

(number  of  dissimilar 
nearest  neighbors) 


3Reference  1. 
Reference  2. 


point,  but  short-range  ordering  persists.    Below  the  transition,  however, 
a  long-range  order  exists  and  the  crystal  exhibits  a  net  magnetization 
over  a  macroscopic  volume  called  a  domain. 

Where  ferromagnetic  ordering  is  increasing  most  rapidly,  we  find  a 
large  specific  heat  and,  as  we  shall  see  later,  a  rapid  length  change. 
The  existence  of  short-range  order  immediately  above  the  transition  is 
shown  dramatically  by  the  large  specific  heat  in  this  region.    For  tran- 
sitions which  occur  at  low  temperatures  where  nonmagnetic  effects  are 

negligible,  the  effects  of  ordering  are  observed  to  about  10T  .    At  tem- 

c 

peratures  much  less  than  T  ,  on  the  other  hand,  the  few  deviations  from 

c 

perfect  ordering  may  be  described  by  "spin  waves"  which  give  rise  to  a 
small  specific  heat. 

Observations  On  Insulating  Ferromagnets 

As  the  experiment  to  be  described  in  Chapter  IV  concerns  an  observa- 
tion on  a  ferromagnetic  salt,  it  is  useful  to  briefly  review  the  available 
information  relating  to  similar  materials. 

There  are  but  a  few  known  ferromagnetic  insulators.     In  most  cases  com- 
parison with  the  theory  is  difficult.    For  example,  gadolinium  trichloride, 
with  a  Curie  temperature  of  2.2  K,  has  large  dipolar  interactions  in 
addition  to  exchange.^    The  same  is  true  for  dysprosium  ethyl  sulphate 
which  has  a  Curie  temperature  of  0.13  K,  so  the  exchange  effects  are 
masked.      Several  workers  have  studied  chromium  tribroraide^'    but  the 
specific  heat  anomaly  at  37  K  is  difficult  to  interpret  because  of  the 
large  lattice  contribution. 

Detailed  studies  have  been  made  of  EuO  (Tc=  69  K)  and  EuS  (T    =  16  K) . 
Both  have  cubic,  nearly  ideal  Heisenberg  lattices  and  the  Eu++  ion  has  spin 


1/2.    The  susceptibility^  '  and  magnetization*^' ^  have  been  measured. 

12 

In  addition,  for  EuO,  the  low- temperature  anisotropy  field,  specific 

heat,**''  ^  thermal  expansion  coefficient,*^  and  induced  magnetostric- 

18  19 
tion     are  known.    Callen  and  Callen      have  applied  the  two-particle 

cluster  approximation  to  EuO  and  EuS.     Isotropic  Heisenberg  correlations 

are  assumed  for  nearest  and  second-nearest  neighbors  only.    The  theory 

is  successful  over  a  wide  temperature  range  on  both  sides  of  the  tran- 

sition. 

Salts  with  the  general  formula  Cu(M+)   (X-), -2H  0,  where  M+  is 
+         +  +  -  -  c        *  Z 

K  ,  NH^  ,  or  Rb    and  X    is  CI    or  Br  ,  are  ferromagnetic  with  Curie 

20 

points  at  liquid  helium  temperatures.    Miedema,  et  al. ,      have  studied 

PI  ?? 

CuCNH^Cl^H^1'^  (Tc  =  0.70  K)  ,  CuRb  C1^2H  0  (T    =  1.02  K) ,  and 

CuCNH^Br^^O25  (T,  =  1.73  K)  as  well  as  CuK  Cl^-gH  0  (T    =  0.88  K)  . 

They  find  the  thermal  properties  of  all  salts  are  identical  apart  from 

a  constant  factor  in  the  temperature.     These  salts  crystallize  in  the 

Ik    ?k  2"5 

tetragonal  system  with  the  space  group  D^h  .     *  r    The  £  axis  is  about  6% 
longer  than  the  a  axis  (a  axis  =  7.45  K  in  CuKgCl^HgO)  so  the  structure 
closely  resembles  body-centered  cubic.    Each  copper  ion  has  eight  nearest 
and  six  next-nearest  copper  ions.    The  unit    cell  is  diagrammed  in  Figure 
1.    There  are  two  molecules  in  the  unit  cell  with  the  copper  ions  at 
(0,0,0)  and  (1/2,1/2,1/2)  positions.    A  water  molecule  is  located  on  each 
side  of  the  copper  ion  in  the  c  direction.    The  four  chlorine  ions  form 
a  rhombohedron  in  the  (001)  plane  about  the  copper  ion.    At  room  temper- 
ature, the  two  magnetic  sites  are  inequivalent  in  the  chloride  salts, 2^ 
but  there  is  only  one  type  of  site  in  Cu(NHk)2Brk«2H  O.2^    For  the  former 
salts,  the  long  diagonal  of  the  chlorine  rhomboid  at  the  cube  center  is 
perpendicular  to  the  long  diagonal  at  the  cube  edge. 

The  salt  CuKgCl^HgO  has  been  the  subject  of  much  study.    Ford  and 


O  x- 

®  H20 
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->  a 


Figure  1.     The  Unit  Cell  of  Salts  Having  the  Formula  CuMgX^HgO 


27 

Jeffries      find  that,  upon  cooling  to  liquid  nitrogen  temperatures, 

perfect  crystals  undergo  a  transition  which  leaves  all  magnetic  sites 

equivalent.     They  suggest  the  long  diagonals  of  the  chlorine  rhomboids 

become  parallel.    Measurements  by  G.  Seidel  and  S.  H.  Choh  reported  by 
27 

Ford  and  Jeffries      suggest  that  the  magnetic  sites  become  equivalent 
because  the  rhomboid  changes  to  a  square. 

Paramagnetic  resonance  absorption  studies22'2^'2^1  show  that  the 

g-value  is  axially  symmetric  about  the  long  diagonal  of  the  chlorine 

32 

rhombohedron.     Paramagnetic  relaxation      gives  approximate  agreement 
with  the  exchange  integral  found  by  the  paramagnetic  resonance  measure- 
ments.    The  susceptibility  is  found  to  obey  the  Curie-Weiss  law  down  to 
32  21 

about  1.6  K      or  1.8  K.        The  zero-field  susceptibility  reaches  a 

value  equal  to  the  inverse  of  the  demagnetizing  factor  in  both  the  c 

and  the  a  directions.    The  magnetization  curves  are  also  similar  in 

both  directions.    A  proton  nuclear  magnetic  resonance  study2^  shows  the 

molecular-field  model  yields  accurate  predictions  for  the  ratio  of  the 

magnetization  to  its  saturation  value.     In  addition,  the  anisotropy 

field  is  found  to  be  quite  small. 

20  21 

Miedema,  et  al. ,  '  find  the  lattice  and  nuclear  hyperfine  con- 
tributions to  the  specific  heat  are  negligible  in  the  range  of  measure- 
ments.   The  specific  heat  is  given  by 

C/R  =  0.262/T2  T  >  I.5  K 

{-O.kSj  log  e  +  0.43  10"^  <  e  <  10_1  T  <  T 

c 

-0.^7  log  £  -  0.05  10"5  <  €  <  lo"1  T  >  T 

c 

where  e  is  the  reduced  temperature  defined  by  e  =  It  -  T  l/T  .  In 

c  1    c " 

addition,  the  specific  heat  of  CuCNH^Br^-  2^0  is  described  by  spin 


C/R 


7 

wave  theory  below  T  =  0.5  T^.    This  salt  was  chosen  for  low- temperature 
analysis  as  the  Curie  temperature  is  relatively  high,  but  the  same 
results  are  expected  to  hold  for  CuKgCl^-2HgO. 

Wood  and  Dalton'^  have  examined  several  critical  properties  of 
CuK2Cl^*2H20  assuming  a  bcc  lattice  with  first  and  second-nearest  neigh- 
bor Heisenberg  interactions  characterized  by  exchange  parameters  and 
J2>  respectively.    It  is  concluded  that  Jg  is  definitely  positive  (ferro' 
magnetic  interaction)  and  J  /J^O^Q.j.     The  high- temperature  series 
yield  values  for  some  critical  constants  shown  in  Table  2.  Non-inter- 
acting spin-wave  theory  is  used  to  reproduce  fairly  accurately  the 
magnetization  and  specific  heat  below  O.k  Tfi  using  estimates  of  and 
J2  obtained  at  higher  temperatures.     It  is  concluded  that  the  exchange 
parameters  are  given  by 

J1/k  =  0.28  +  0.01  K 

Jg/k  =  0.056  +  0.003  K. 

To  summarize,  we  note  that  the  four  salts  given  above  with  the 

general  formula  CuMgX^'2HgO  give  useful  information  for  the  following 
20,27,33 

reasons: 

1.  The  magnetic  interactions  are  due  largely  to  first  and  second- 
nearest  neighbor  Heisenberg  exchange. 

2.  The  Curie  temperatures  are  low  enough  that  the  lattice  contribu- 
tion to  the  specific  heat  and  the  expansion  coefficient  is  negligible. 

3.  The  copper  ion  has  spin  1/2,  so  high- temperature  expansions  for 
the  spin  1/2  model  can  be  used. 

k.    The  lattice  is  uncomplicated,  and,  at  low  temperatures,  all 
magnetic  sites  are  equivalent. 

Subsequent  chapters  will  consider  the  specific  heat  and  the  thermal 


Table  2.    Comparison  of  Experiment  and  Theory  for  Some 

Critical  Constants  of  CuKjCI.  *2H  0 

2    k  2 


3 

Property  Experiment  Theory 

T  0.88  0.85 
c 

Curie-Weiss  constant  1.2  +  0.1b  1.29 

Energy  difference  between 
paramagnetic  and  ferro- 
magnetic states  at  T  =  0  5.50°  5.3^ 

C^/R,        T   >  2Tc  2.62C,  2.56d  2.36 

(Sm  -  S  )/k  0.222C  0.229 

(E    -  E  )/k  0.384c  0.419 

w  C 


a 

Reference  33- 

^Reference  21. 

°Reference  20. 
d 

Reference  32. 


expansion  coefficient  near  the  critical  temperature.    Unfortunately,  the 
critical  behavior  of  these  functions  is  not  easily  determined  experimen- 
tally.   Many  researchers  note  difficulty  in  choosing  one  form  over  an- 
other.    The  situation  is  often  further  clouded  by  corrections  for  lattice 
effects,  for  the  shift  in  exchange  parameters  with  volume,  and  for  round- 
ing of  the  lambda  peak. 

20  21 

As  noted  previously,  Miedema,  et  al. .     '      have  found  that  the  lattice 
contribution  to  the  specific  heat  in  CuKgCl^HgO  is  negligible  over  the 
range  of  their  measurements  (below  k.2  K) .    To  get  an  order-of-magnitude 
estimate  of  the  lattice  contribution  to  the  thermal  expansion  coefficient, 
0L,  assume  a  Gruneisen  behavior  for  the  lattice  so  that 


where 


3V 


y  -  .  <*  I2S  9 
'l  d  log  V 


and  where  e  is  the  Debye  temperature.    Assume  v^~>~l  as  in  the  alkali 

halides  for  which  the  expansion  coefficients  have  been  measured  at  low 

temperatures.        Assume  also  that  a  similar  relation  holds  for  the 

magnetic  contribution  to  the  thermal  expansion  coefficient,  a  ,  as  in 

m 

17 

the  model  of  Argyle,  et  al. .      described  in  Chapter  II.    Assume  also 
that 

y  .    .    d    lSB    J       =  10/, 

'm  d  log  V  iU° 

as  suggested  by  Bloch      for  systems  coupled  by  superexchange  inter- 
actions. Then 


It  is  thus  reasonable  to  expect  that  the  lattice  contribution  is  negli- 
gible in  the  thermal  expansion  coefficient  as  well  as  in  the  specific 
heat  of  CuK  Cl^'ay)  below  k.2  K. 

The  shift  of  the  exchange  parameters  with  volume  may  be  most  simply 
expressed  as  a  shift  of  the  effective  Curie  point  with  temperature  as 
follows: 

dT 

AT     ■    -rjr-  AV 
c  dV 

d  log  T 

=    T     (3  —  AT, 

c       d  log  V 

where  £  is  the  volume  thermal  expansion  coefficient.    Heller  and 
56 

Benedek       calculate  a  temperature-dependent  Neel  temperature  for  MnF^. 

An  alternate  approach  to  the  correction  for  the  shift  of  the  ex- 
change parameters  with  volume  assumes  a  constant  critical  temperature  but 
shifts  the  temperature  scale  slightly.    Assume  for  definiteness  a 

negative      _  *°8  Tc    and  a  positive  expansion  coefficient.    Then  as  T 
d  log  V  C 
is  approached  from  below,  the  sample  expands  and  T£  decreases.    For  any 

temperature,  however,  the  distance  to  the  effective  transition  tempera- 
ture is  greater  than  |t  -  T  | .    The  same  result  holds  for  T>  Tc.  Hence 
Yamamoto,  et  al. ,      make  the  simplest  correction  and  replace  |T  -  TcJ 
with  |T  -  T  |    +  A. 

To  get  an  approximate  value  for  the  change  of  Tc  in  CuKgCl^'Sa^O, 

set    d  lo8  Tc  =-io/3  as  before  and  use  0-\^3a££\^  1.5x10  5K  *|log  c|  near 
d  log  V 

T  . 
c 

Then  ATC^0.2    uK  when      e  =  lo"5 

c  r 

and  AT  C±*>  8    uK  when      €  =  10*1 . 

c 


11 

A  shift  of  this  magnitude  is  not  significant  as  we  know  Tc  only  to  1  mK. 

Cadieu  and  Douglass^  define  the  supercritical  region  as  the  temper- 
ature range  near  T    for  which  the  temperature-dependent  coherence  length 

c 

is  comparable  with  the  spin-fluctuation  mean  free  path.    The  supercritical 

behavior  begins  at  a  temperature  for  which  the  correlation  length 

describing  the  spin  fluctuation  spectrum  becomes  comparable  with  the 

mean  free  path  corresponding  to  some  process  which  tends  to  destroy 

ordering.    Let  €    be  the  reduced  temperature  giving  the  onset  of  the 
s 

supercritical  region.    Cadieu  and  Douglass  find  for  Gd,  with  Tc  =  292  K, 

-0  72  -1 
€  =    constant  x  T  , 

a 

where  T  is  an  indicator  of  impurity.    Teaney*-5  has  attributed  the  severe 
rounding  of  the  specific  heat  anomaly  of  EuO  to  the  presence  of  filamen- 
tary inclusions  of  tantalum  metal. 
^9  kO 

Domb    *      has  undertaken  a  comparison  of  graphs  which  contribute  to 

the  partition  function  of  an  infinite  assembly  and  a  finite  assembly 

wrapped  on  a  torus.    This  leads  to  an  estimate  of  the  specific  heat  of 

finite  two-  and  three-dimensional  Ising  models.    The  reduced  temperature 

at  which  the  specific  heat  reaches  a  finite  maximum  is  given  by  M 

where  the  cubic  assembly  has        spins.    An  order-of -magnitude  estimate 

of  €    is  therefore  given  by 
s 

8  » 

where  N  is  the  number  of  magnetic  ions  in  the  assembly. 

After  corrections  for  the  lattice  effects  and  for  the  volume  shift 
of  the  Curie  point  have  been  made,  a  form  is  sought  to  describe  the 
specific  heat  or  expansion  coefficient  outside  the  supercritical  region 
(if  one  is  observed) . 


12 

Teaney*^  reports  that  preliminary  analysis  of  the  specific  heat  of 
EuO  suggests  a  logarithmic  dependence  so  that 

C±  =  A+  log  €  +  B+, 

where  (+)  and  (-)  refer  to  the  critical  region  above  and  below  the 

transition  respectively. 

37 

Yamamoto,  et  al. .      have  analyzed  a  number  of  antiferromagnets  as 
well  as  ferromagnetic  Sm,  Tb,  Gd,  and  p-UH^  in  terms  of  a  "pinched 
logarithmic"  function  given  by 

C+  -  A  log  (e  +  A)  +  B+. 

Unfortunately,  the  possibility  of  some  other  kind  of  singularity  remains 

in  cases  other  than  CoClg'SHgO  and  Tb. 

20 

Miedema,  et  al. ,      studied  the  specific  heat  of  CuKgCl^^gO  in  de- 
tail between  0.9  Tc  and  1.1  Tfi.    As  pointed  out  above,  they  find 

C+  «  A  log  e  +  B+. 

This  is  consistent  with  the  "pinched  logarithmic"  form  since  the 
correction  for  the  volume  dependence  of  T    is  quite  small. 

Chapter  II  indicates  how  linear  thermal  expansion  coefficient  data 
can  shed  light  on  certain  properties  of  the  ferromagnetic  insulator. 
First,  it  is  shown  that  the  temperature  derivative  of  certain  spin 
'correlation  functions  describes  the  expansion  coefficient  and  the  specific 
heat.     If  the  specific  heat  is  proportional  to  the  expansion  coefficient, 
further  information  may  be  obtained  about  these  spin  correlation  func- 
tions and  the  stress  dependence  of  the  transition  temperature.  Next, 
thermodynamics  of  anisotropic  materials  near  a  lambda  point  suggests 
another  method  to  determine  the  stress  dependence  of  T  .    Finally,  the 
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expected  behavior  of  the  thermal  expansion  coefficient  and  the  specific 
heat  is  discussed. 

Chapter  III  describes  the  experimental  procedure  and  the  apparatus 
we  have  used  to  measure  the  linear  thermal  expansion  of  a  small  crystal 
of  ferromagnetic  CuK2Cl^-2H20  at  liquid  helium  temperatures. 

Chapter  IV  gives  the  method  of  analysis  and  the  results  of  our 
thermal  expansion  coefficient  measurement.     It  is  concluded  that  both 
significant  spin  correlation  functions  have  the  same  temperature  de- 
pendence well  above  the  transition.    Furthermore,  the  stress  dependence 
of  the  transition  temperature  found  in  two  ways  agrees  within  experi- 
mental error. 


CHAPTER  II 
THEORY 

The  Magnetic  State 

The  magnetic  phenomena  in  solids  are  cooperative  effects  which  are 
manifestations  of  correlations  between  magnetic  ions.    A  series  of  papers 
by  P.  W.  Ander  son  has  shed  considerable  information  on  the  nature  of 

these  correlations,  and  the  brief  discussion  which  follows  is  taken  from 
this  work. 

Attempts  to  describe  the  magnetic  state  in  insulating  solids  lead 
to  a  paradox.    Whereas  interactions  between  magnetic  ions  are  essential 
to  cooperative  behavior,  the  electrons  must  be  localized  as  conduction 
is  not  observed.    Suppose,  for  example,  the  magnetic  ion  is  Cu**  with  a 
single  d  hole.    Consider  a  single  d  hole  allowed  to  wander  throughout  the 
lattice.    We  may  put  it  at  the  bottom  of  the  d  band.    When  the  hole  is 
in  the  vicinity  of  another  ion,  it  is  repelled  by  the  Coulomb  energy  of 
the  localized  d  hole.     In  dilute  materials  for  which  the  energy  bands 
are  narrow,  the  decrease  in  kinetic  energy  is  outweighed  by  the  Coulomb 
correlation  energy,  U,  so  the  hole  remains  localized. 

If  a  solid  were  simply  a  lattice  of  non-overlapping  free-ion  wave- 
functions,  paramagnetic  behavior  would  be  observed  at  all  temperatures. 
The  interactions  which  lead  to  spin  correlations  are  of  electrostatic 
origin,  however,  and  are  a  function  of  the  wavef unctions  of  the  magnetic 
ions  and  the  non-magnetic  ligands  in  the  solid.    Anderson  has  described  a 
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perturbation  procedure  which  allows  one  to  discuss  the  magnetic  effects 
without  becoming  involved  in  the  complex  problems  of  wavefunctions  in 
solids . 

What  is  desired  is  a  procedure  for  separating  the  problem  into  two 
parts  such  that  the  part  containing  the  magnetic  interaction  is  not  a 
large  perturbation  on  the  other  part.    The  He i tier-London  treatment  re- 
quires that  the  interaction  of  magnetic  ions  with  the  ligands  be  treated 
as  part  of  the  perturbation.     This  is  a  much  larger  perturbation  than 
the  magnetic  effects  alone. 

Assume  only  one  d  hole  per  magnetic  site  for  def initeness .  Suppose, 
with  Anderson,  we  consider  the  diamagnetic  lattice  formed  by  filling  the 
d  hole  in  every  magnetic  ion  but  one.    Adjust  the  nuclear  potential  so 
the  single  hole  sees  the  correct  potential  as  it  moves  through  the  lattice. 
This  is  the  ligand  field  theory  aspect  of  the  problem  and  involves  the 
complete  solution  of  the  lattice  exclusive  of  exchange  effects.     The  one- 
hole  wavefunction  may  be  a  Bloch  running  wave  orthogonalized  to  the 
magnetic  ion  core,  but  since  the  ligand  engages  in    covalent  bonds  with 
the  magnetic  ion,  the  core  wavefunction  contains  a  mixture  of  d  functions. 
The  Bloch  wave  must  then  contain  an  antibonding  admixture  of  the  ligand 
function.    This  is  the  mechanism  which  gives  rise  to  superexchange  via 
non-magnetic  ions,  although  this  treatment  shows  that  once  the  non-magnetic 
problem  of  the  lattice  without  exchange  effects  is  solved,  the  ligands 
need  not  be  considered  further. 

Anderson  has  shown  that  the  effect  of  the  ligands  can  be  summarized 
in  transfer  integrals  b    ,  (R  -  R')  which  allow  transitions  between 
magnetic  ion  sites  R  and  R1  and  states  n  and  n'  without  change  of  spin. 
All  other  effects  follow  from  many-electron  interactions. 
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The  b  integrals  allow  a  hole  to  decrease  its  kinetic  energy  by  making 

virtual  transitions  to  other  sites.    This  decrease  in  kinetic  energy  is 

2 

given  roughly  by  2b  /U,  where  U  is  the  Coulomb  correlation  energy  discussed 
earlier.    Clearly,  the  transition  is  only  allowed  if  the  spins  at  the  two 
sites  are  opposite,  since  the  b's  do  not  cause  a  change  in  spin.  Hence 
we  can  represent  the  energy  shift  by 

I*b2/U(-1A  +  S-S'). 

This  is  not  a  true  exchange  effect  as  it  occurs  for  antiparallel  spins, 
although  it  is  useful  to  ignore  this  distinction  and  label  the  effect 
kinetic  or  antiferromagnetic  superexchange . 

Other  spin-dependent  inter-ion  effects  follow  using  the  Hartree-Fock 
machinery  with  the  wavef unctions  obtained  from  the  solution  of  the  ligand 
field  aspect  of  the  problem.     The  most  significant  of  these  effects  is  the 
familiar  ferromagnetic  potential  exchange  given  by 

•^.(R  -  R')(iA  +  s-s'), 

where  J  is  the  Coulomb  self -energy  of  an  "overlap  charge  distribution" 
using  states  n  and  n'. 

Although  other  effects  can  be  treated  with  the  present  method,  anti- 
ferromagnetic  kinetic  exchange  resulting  from  virtual  transfer  between 
magnetic  ions  and  ferromagnetic    Heisenberg  potential  exchange  usually 
dominate  in  materials  for  which  magnetic  effects  are  due  to  a  single  d 
hole  or  electron. 

One  concludes  that  the  effective  Hamiltonian,  H,  for  insulators  with 
spins  coupled  by  superexchange  interaction  is 

H  =    Z  J(R  -  R«)  s-S' 

magnetic  ionv  ' 

sites  R,  R' 
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where  constant  terms  are  neglected  and  the  sign  of  the  exchange  parameters 
Is  positive  or  negative  as  kinetic  or  potential  exchange  dominates.  This 
is  the  familiar  Heisenberg  Hamiltonian. 

Spin  Correlation  Functions 

The  treatment  of  Anderson  justifies  the  Heisenberg  Hamiltonian  for 
systems  coupled  by  superexchange  interactions.    However,  an  efficacious 
procedure  for  determining  the  thermal  expansion  coefficient  from  this 
Hamiltonian  must  entail  further  assumptions  as  there  is  no  exact  solution 
of  a  three-dimensional  Heisenberg  system.    A  useful  statement  of  the  prob- 
lem has  been  offered  by  Callen  and  Callen^  who  have  given  expressions 
relating  the  macroscopic  strains  to  products  of  microscopic  coupling  con- 
stants and  spin  correlation  functions.     Their  development  includes  single- 
ion  crystal  field  effects  and  two-ion  interactions.     The  strains  are 
characterized  entirely  by  three  spin  correlation  functions:     the  one-ion 
longitudinal  correlation  function  which  varies  as  <(S*)S>  (subscripts  label 
magnetic  ions),  the  two-ion  isotropic  correlation  function  which  varies  as 
<S£-Sg>»  and  the  two-ion  longitudinal  correlation  function  which  varies  as 

<isK  -  jvv>- 

It  is  useful  to  apply  the  procedure  described  by  Callen  and  Callen  to 

the  tetragonal  system  appropriate  for  CuKgCl^HgO.    We  assume  there  is  no 

pi 

easy  or  preferred  direction  of  magnetization,      so  the  magnetization  with- 
in a  domain  is  equally  likely  to  lie  in  any  direction.    Furthermore,  the 
one-ion  longitudinal  correlation  function  is  constant  for  spin-1/2  systems 
and  therefore  does  not  contribute  to  the  expansion  coefficient.    Hence  the 
linear  expansion  coefficient,  a,  measured  in  a  direction  specified  by  the 
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unit  vector  B  is 

z  [k'bb  +kJ   /n      _  v„  _  a  <s2sz  -  -s  -s  >, 

f  ,g      f*  *  z       f8  <BX  +  VBZ3  §! ■    «  8      3  f    8  * 

where  the  K's  are  products  of  elastic  constants  with  two-ion  magneto- 
elastic  coupling  coefficients.  For  a  along  the  a  axis,  only  one  com- 
ponent of  B  is  nonzero.    We  have  then  that 

a  =  Z  K.    |r  <S  -S  >. 

Since  all  magnetic  sites  are  equivalent  and  first-  and  second-nearest 
neighbors  dominate  the  interaction,  the  sum  has  only  two  terms: 

am"zihh  ^i'V  +  Nz2K2    <si-s3>*  U] 

where  N  is  the  number  of  ions  per  unit  volume,        is  the  number  of 
nearest  neighbors  (z^  =8),  zg  is  the  number  of  next-nearest  neighbors 
(z2  =  ^) »  anc*  *ons  two  and  three  are  representative  nearest  and  next- 
nearest  neighbors  of  ion  one. 

Since  the  K's  are  functions  of  the  stresses  only,  Eq.  [l]  is  deriv- 
able from  a  Gibbs  free  energy  given  by 

G  =  J^YT)  +  J2g2(J2/T), 

with 

1  *1 
*1  "      V    dcx  » 

1  ^2 

and      K2  =  -7  sr* 

and  where  o  is  the  component  of  stress  in  the  a  direction,  V  is  the  molar 
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volume,  and  the  exchange  parameters  are  functions  of  the  stresses  only. 
Then  finally  we  have 

rt  =    .   1   1    <L     <rs  -S  >  -   -  — -   —    <S1       >*  [21 

It  is  interesting  that  the  magnetic  contribution  to  the  specific 
heat  at  constant  stress,  C^,  is  given  by 

Suppose  that  the  a  axis  linear  expansion  coefficient  and  the  constant- 
stress  specific  heat  are  found  to  be  proportional.    Let  P  =  V  Ct/C^.  Then 
one  or  both  of  the  following  is  true: 

d  In  J         d  In  J 

1         P  -  1  -   — 

do  dc 

b  [ln(z1  +  rz2)  +  In  J  ^         bin  (z^  +  z^) 

«  _  =  §5 

where  r  =  Jg^i  (about  °*3  for  CuKgClr *2HgO),"  We  assume  the  transition 

temperature  scales  with  the  interaction  energy.  Then  we  have 
d  In  T 

*  -  — sr*  •  t*i 

2.    The  two  spin  correlation  functions  have  proportional  temperature 
derivatives.     In  this  case,  let 

h  <si's3>  ■  M  h  <YS2>' 

where  M  is  a  constant.  Then 


20 

d  ln(z1J1  +  Mz2J2) 
p  =  gi  • 

Callen  and  Callen  have  applied  cluster  approximation  methods  to  find  that 

the  correlation  functions  have  approximately  equal  temperature  derivatives 
19 

above  T^/^.        These  calculations  are  not  expected  to  apply  in  the 
critical  region,  but  they  suggest  that  well  above  the  transition  M  =  1 
so  that 

p       d  In  T 

P  =   c  . 


Thermodynamics 

It  is  necessary  to  prelude  the  study  of  the  thermodynamics  of  the 
lambda  transition  with  a  discussion  of  the  critical  behavior  of  the  com- 
pressibility.   The  model  we  use  is  not  adequate  for  CuKgCl^HgO,  but 
order-of-magnitude  results  are  sufficient  to  show  that  the  compressibility 

is  constant.    Then  the  thermodynamics  of  anisotropic  media  near  a  critical 

dT 

point  can  be  applied  to  suggest  another  method  to  determine  *r~  .  A 
formula  will  also  be  derived  for  the  shift  of  the  Curie  point  with  strain. 

1  7 

Argyle,  Miyata,  and  Schultz      have  considered  a  model  which  assumes 
the  Hamiltonian  is  the  sum  of  a  Gruneisen  lattice  terra  plus  a  magnetic 
term  involving  exchange  constants  which  all  exhibit  the  same  volume  de- 
pendence.   The  differentiations  to  find  the  compressibility,  K,  give  a 
lattice  term,  K  ,  and  a  magnetic  compressibility,  K  ,  related  by 

K=  ^ra'^VV1    =K1<1  +  W'1. 

where 
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K  -1  =  -T/VC 
m  m 


f  \ 
d  log  J 
d  log  V 


2 

+  term  negligible  near  T  .  [6] 


c 


We  estimate  the  ratio  K  /K  .    Assume  j  is  approximately  -10/3 

1    m  a  log  v 
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as  suggested  by  Bloch      for  systems  coupled  by  superexchange  interactions. 

!■=  —  I46 

Reviews  by  Huntington  J  and  Aleksandrov  and  Ryzhova      have  attempted  to 

collect  all  known  data  on  the  elastic  constants  of  crystals.    The  latter 
review  is  complete  through  I960.     Information  about  the  compressibility 
of  crystals  similar  to  CuK  Clk'2H  0  is  seriously  lacking.    Fjrom  the  com- 
pilation of  Huntington,  we  take  a  typical  value  for  the  alum  salts.  Hence 

-12  2 

we  assume  K^/~\J  7x10      cm  /dyne.     In  the  critical  region  where  de- 
creases most  rapidly,  the  asymptotic  part  of  the  magnetic  contribution 
to  the  specific  heat  is      0.l467R|log  e|.     In  this  crude  approximation, 
we  use  Eq.  [6]     to  find 

■50000. 


and 


K./K     =  -  1  when  e'^lO 

1  m 


_k  -k 

K./K     =  -  10  when  €^X/10  . 

1  m 


This  result  has  two  consequences. 

1.  The  compressibility  may  be  considered  constant  and  equal  to 
that  measured  away  from  the  critical  region. 

2.  A  more  fundamental  consequence  concerns  the  nature  of  the  trans- 
ition, for  when  Ik./K  I  >  1,  the  total  compressibility  is  negative.  This 

■  1  m 

k7 

leads  to  a  first-order  transition.        The  finite  size  of  the  magnetic 

domain  will  limit  the  divergence  of  the  specific  heat  even  in  perfect 

crystals  when  the  reduced  temperature  is  on  the  order  of  10  ^ ^  Hence, 

the  lambda  nature  of  the  transition  is  preserved. 
1+S 

Janovec      has  extended  to  anisotropic  dielectrics  the  Buckingham  and 
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Fairbank      treatment  of  the  thermodynamics  of  the  lambda  transition.  To 
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proceed  with  an  analogous  development  for  f erromagne ts ,      a  relation 
among  partial  derivatives  is  needed. 


Suppose  W  is  a  function  of  seven  of  the  eight  variables  y,z,  and  x^ 
where  i  =  1, ...,6.    Let  us  adopt  the  following  conventions: 

1.  Use  x  to  represent  the  set  of  six  variables  x^  i  =  i, ...,6. 

2.  Use  Xj^  to  represent  the  set  of  five  variables  x^,  j  =  1,...,6, 
except  j  £  i. 

Then 


dW  =  ,Z 
k 


K 

I  / 


dx,  + 
k 


/ 

dy 


dy 


or 


dW  = 


/  xj/i> 


dy  + 


2w 

dz 


Xj/i'y 


dz  +  Z 
k^i 


3w 


dx. 


Considering 


xi  =  xi<xj/i>y>z>> 


dx.  = 


k^i 


V 


dx,  + 
k 


Cj/i,k'y>z 


\  / 


dy  + 


Xj/i,Z 


dz. 


'j/i': 


Equating  the  value  of  dW  and  setting  d    =  o  and  dx.,.  =  o,  we  find  the 

z  j  /  i 


following  identity: 


dW 
dy 


/ 


X 


d^ 
dy 

\  / 


x.  ,z 


x . 
i 


x.  , .  ,y 


dy 


[7] 


xj/i'z 


It  is  useful  to  introduce  the  "neighborhood  temperature,"  t,  given  by 


T  -  T 


The  properties  of  t  are  discussed  by  Buckingham  and  Fairbank 


We  will  consider  Tc  a  function  of  the  stresses  only. 
The  specific  heat  at  zero  stress  is 


C    =  T 
a 


as 

dT 


\  / 

Using  Eq.  [7],  we  find 


dT 

\  I 


I 

as 
ai 
\  1 


ffj/i'c 


as 
5, 


/ 


cj/i'T 


ai 


so  that  we  have  for  all  temperatures 


C    =  T 
a 


as 
aT 

c 

\  I 


TV  a. 


ai 

 c 

aa. 


'j/i 


[8] 


where 


°i  - 


i 

aT 


and  eA  is  the  strain  function. 

Assuming  the  shape  of  the  transition  is  not  a  strong  function  of 
stress,  the  rate  of  change  in  entropy  with  respect  to  stress  at  a  distance 
t  from  the  transition  is  small.    Since  T,  is  a  function  of  stress  only, 
neither  is  the  rate  of  entropy  change  with  respect  to  T.  large.  Hence 

1  as 

T  I  aic  j   is  sma11  and  constant  compared  with  the  large  and  rapidly  vary- 
ing expansion  coefficient  and  specific  heat.    Thus,  to  a  good  approxima- 


C    =  constant  -  T  Va 
a  c  i 


ac 


'j/i 


near  T  . 
c 


[9] 
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The  constant  of  proportionality  between  the  asymptotic  parts  of  the 
specific  heat  and  the  linear  expansion  coefficient  gives  the  stress  de- 
pendence of  the  transition  temperature. 

The  dependence  of  the  transition  temperature  on  the  strains  is  of 
interest.    We  use  Eq.  [7]  again  to  find 


1  \ 

dT 

dT 

a 

i  I 

dT 

\  / 

°j/i'T 

1  i 

or 


ai  = 


dT 

c 

ST 


-1 


-  s 


ii 


3T 


-1 


Si  / 


where  the  s^  are  the  moduli  of  compliance  for  the  crystal  given  by 


'ki 


°j/i'T 


Let  us  estimate  the  value  of  each  term  in  this  relation.    First,  we  find 


dT 

c 

s; 

\  / 


-1 


5T 


d  log  T 


d  log  V 


for  isotropic  materials.  Assume 
d  log  V    3 


as  suggested  by  Bloch 
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c 

Si 

\  / 


-1 


Then  we  have 
-1 


^>   -  0.1  K 
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Take  as  an  order-of-magnitude  estimate  for  s^  the  room  temperature  com- 
pressibility  of- alum  '  divided  by  three.    Our  measurements  indicate 


dTc  1 

I  =  -  1.5x10      K  cni  /dyne 

1 


for  the  a  direction,  so  we  estimate 

■1 


c 


ii 


do, 


-0.2K 


'j/i 


■6„-l, 


The  asymptotic  value  of  a  in  the  a  direction  is  about  5x10  K  |  log  c|, 
so  over  the  range  of  measurements, 


a   <  3xlo"5K_1. 


Hence       may  be  neglected  with  respect  to  the  other  terms  and  we  have 


t 

c 


ii 


c 


[10] 


'j/i 


Etl«  [10]  gives  the  dependence  of  the  interaction  energy  (proportional 
to  Tc)  on  the  lattice  parameters.    These  quantities  can  be  computed 
from    the  less  interesting  stress  dependences  of  the  transition  tempera- 
ture only  if  the  moduli  of  compliance  are  known.     It  is  unfortunate  that 
these  are  not  available  for  any  salts  having  the  general  formula 
CuM2X4-  2H20. 

Predicted  Behavior  of  the  Thermal  Expansion  Coefficient 

It  is  desirable  to  compare  the  measured  expansion  coefficient  with 
the  behavior  predicted  by  the  theory.     This  comparison  is  possible  for 
CuKgCl^.  2H20  in  three  distinct  temperature  regions  as  follows: 
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1.  The  spin  wave  approximation  has  been  successfully  applied  to 
the  specific  heat  of  CuCNH^Br^  2H2020  at  temperatures  below  0-5  T  . 
Our  data  for  CuKgCl^aigO  cannot  be  compared  with  this  theory  as  measure- 
ments do  not  extend  to  sufficiently  low  temperatures. 

2.  The  magnetic  contribution  to  the  specific  heat  is  expected^'^1 
to  be  proportional  to  T     well  above  the  Curie  point.    This  prediction 
is  verified  by  Miedema,  et  al..20  and  Van  den  Broek,  et_al . 52    The  argu- 
ment of  Callen  and  Callen  discussed  previously  suggests  that  the  expan- 
sion  coefficient  will  vary  as  the  specific  heat  in  this  temperature  range. 

3.  In  the  critical  region  (reduced  temperatures  less  than  about 
0.1),  no  form  for  the  thermal  expansion  coefficient  is  predicted  unam- 
biguously by  the  theory. 

Domb  has  recently  published  a  review  of  critical  behavior  of  magnetic 
52 

systems.        Numerous  references  to  previous  work  appear;  they  will  not 
be  repeated  here.    Several  general  considerations  discussed  by  Domb 
should  be  mentioned. 

1.  The  critical  behavior  of  the  thermodynamic  functions  is  more 
strongly  influenced  by  the  dimensionality  of  the  system  than  by  details 
of  the  lattice  structure.    For  example,  the  Curie  point  of  a  two-dimen- 
sional Heisenberg  ferromagnet  is  at  absolute  zero,  while  T^"\-AJ/k  for  all 

c 

three-dimensional  systems. 

2.  The  behavior  of  the  susceptibility  above  the  transition  is  pre- 
dicted theoretically  and  there  is  reasonable  experimental  verification. 
The  specific  heat  is  much  more  difficult  to  predict  and  depends  more 
strongly  on  the  model  chosen. 

3.  The  divergence  of  the  specific  heat  above  the  transition  is 
somewhat  slower  for  Heisenberg  systems  than  for  Ising. 


Domb  suggests  that  for  a  Heisenberg  ferromagnet  with  spin  1/2,  the 
specific  heat  diverges  slowly  and  may  be  logarithmic  for  all  lattice 
struc tures. 
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Baker,  et_al . ,      have  analyzed  the  specific  heats  of  several  spin- 
1/2  Heisenberg  lattices.    For  the  bcc  lattice  they  find 

C  =  (Tc/T)2[A  -  B(l  -  Tc/T)0-20]      0.05  <    €    <  0.3,  X  >  *  . 

If,  however,  one  allows  the  critical  temperature  as  determined  by  the 
specific  heat  to  differ  from  that  determined  from  the  susceptibility, 
a  logarithmic  function  is  also  reasonable. 


CHAPTER  III 
EXPERIMENTAL  PROCEDURE  AND  APPARATUS 
The  Sample 

The  CuKgCl^'f^O  sample  upon  which  our  thermal  expansion  measurements 
were  made  was  prepared  by  Dr.  George  S.  Dixon,  Jr.  at  the  Oak  Ridge 
National  Laboratory,  who  studied  an  adjacent  section  of  the  same  crystal 
and  found  that  piece  to  be  a  single  crystal  and  of  high  quality.  Our 
piece  has  visible  faults,  however. 

As  the  saw  cut  on  the  material  we  received  lay  in  the  (100)  plane, 
we  chose  to  make  our  measurements  perpendicular  to  this  cut.    After  per- 
forming the  necessary  cuts  with  a  Lastec  saw  using  a  non-abrasive  wire 
moistened  with  water,  the  sample  height  was  only  O.58  cm.     The  two  faces 
in  the  (100)  plane  were  polished  by  moving  the  crystal  about  on  a  very 

slightly  moistened  cloth  stretched  across  a  smooth  surface.     These  faces 

2 

are  irregular  in  shape  and  have  an  area  of  approximately  1  cm  . 

Our  experience  with  hydrated  crystals  cautions  against  evacuating  the 
space  around  such  a  material  at  temperatures  greater  than  0  C.    Care  was 
exercised  to  see  that  the  sample  chamber  was  at  atmospheric  pressure  un- 
less it  was  well  below  the  freezing  point  of  water.    Furthermore,  these 
salts  tend  to  deteriorate  if  exposed  to  air.    We  therefore  stored  our 
crystals  in  vacuum-pump  oil. 
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As  CuK2C1^2H20  reacts  with  the  copper  sample  chamber,  1/2  mil  mylar 
sheets  were  sandwiched  between  the  sample  and  the  copper.    These  sheets 
are  quite  thin  and  do  not  affect  the  measurements.    A  coating  of  varnish, 
grease,  or  heavy  oil  might  apply  stress  at  low  temperatures  and  was  not 
used. 

The  Cryogenics 

The  cryostat  is  similar  to  that  described  by  Walsh.51*    It  has  been 
described  by  Straty,55  although  a  number  of  small  changes  have  been  made. 
Figure  2  shows  the  cryostat  schematically.    Both  dewars  are  sealed  at  the 
top  with  rubber  O-rings  so  the  nitrogen  can  be  pumped  below  its  triple 
point  as  will  be  discussed  later,  and  the  helium  boil-off  can  be  recovered. 
The  tubulation  extending  into  the  liquid  helium  bath  enters  a  vacuum-tight 
cylinder  enclosing  a  container  of  liquid  He4  (the  He4  evaporator)  and  an- 
other evaporator  for  liquid  He3.    The  sample  chamber  is  attached  to  the 
He5  evaporator  and  can  be  cooled  to  nearly  0.3  K. 

The  brass  cylindrical  container  enclosing  the  evaporators  and  the 
sample  chamber  is  attached  to  a  flange  on  the  cryostat  stem  using  a  vacuum- 
tight  gasket  of  0.075  cm  diameter  pure  indium  wire.    After  the  initial 
liquid  helium  transfer,  the  He4  exchange  gas  in  this  container  is  evacuated 
by  pumping  for  several  hours  with  a  Consolidated  Vacuum  Corporation  model 
PMCS-2C  oil  diffusion  pump  backed  by  a  Welch  model  11*00  pump. 

The  He4  evaporator  is  filled  from  the  bath  via  a  modified  Hoke  valve 
operated  by  a  long  shaft  extending  through  the  top  flange  of  the  cryostat. 
The  evaporator  temperature  can  be  maintained  at  approximately  1  K  for 
over  2k  hours  by  pumping  on  the  enclosed  helium  with  a  Welch  model  li«)2 
pump  or  a  high-volume  Kinney  model  KC-k6  pump.    A  radiation  shield  attached 


Figure  2.    Schematic  Drawing  of  the  Cryogenics 


A.  Flexible  stainless  steel  tubing  for  returning  He 

B.  Flexible  stainless  steel  tubing  to  sample  chamber  pump 

C.  Flexible  stainless  steel  tubing  to  manometry  system 

D.  Flexible  stainless  steel  tubing  to  He^  pump 

E.  Quick  vacuum  coupling  for  exchange  gas  pumping  line 

F.  Rubber  hose  for  He    evaporator  pumping  line 

k 

G.  Rubber  hose  to  He    recovery  system 

H.  He^  diffusion  pump 

I.  Rubber  inner  tube 

J.    N    fill-pump  line 

K.    N2  bath 
k 

L.    He  bath 
k 

M.    He  evaporator 

N.     Radiation  shield 
3 

0.    He  evaporator 

P.    Sample  chamber  exchange  gas  container 
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to  the  underside  of  the  He^  evaporator  surrounds  the  He^  evaporator  and 

the  sample  chamber. 

The  He'  refrigeration  system  was  designed  for  continuous  operation. 

3  5 

The  evaporator  holds  approximately  1  cm  of  liquid  He  .  A  National  Re- 
search Corporation  model  B-2  oil  diffusion  pump  backed  by  a  Welch  model 
lk02  pump  with  an  oil  shaft  seal  for  closed  system  operation  can  reduce 
the  evaporator  and  the  sample  chamber  temperature  to  nearly  0.3  K.  The 
return  gas  is  forced  to  pass  through  a  trap  immersed  in  liquid  nitrogen. 
This  trap  consists  of  a  stack  of  discs  cut  from  brass  screen  followed  by 

a  container  of  Linde  13X  molecular  sieve  pellets.    The  return  line  is  con- 

k 

stricted  as  it  passes  by  the  He  evaporator.  This  constriction,  made  by 
feeding  a  10  mil  nichrome  wire  through  about  10  cm  of  12-mil-i.d.  stain- 
less steel  capillary,  reduces  the  pressure  of  the  returning  He'  which  has 

It 

been  liquified  in  the  constriction  through  thermal  contact  with  the  He 
evaporator. 

The  sample  chamber  is  bolted  directly  to  the  bottom  of  the  He'  evapor- 
ator, the  abutting  faces  being  lightly  greased  to  improve  thermal  conduc- 
tion.   A  copper  vacuum-tight  container  encloses  the  sample  chamber.  The 
vacuum  seal  is  again  made  using  a  gasket  of  pure  indium  wire.    The  sample 
chamber  may  be  evacuated  through  a  0. 127~cm-o. d.  stainless  steel  capillary. 
If  all  systems  perform  properly  at  k-  K,  He'  exchange  gas  is  introduced  in- 
to the  sample  chamber. 

Length  changes  of  the  sample  are  measured  using  a  capacitance  tech- 
nique discussed  later.    The  two  capacitor  leads  and  the  two  leads  from  a 
germanium  resistance  thermometer  are  coaxially  shielded.    There  is  a 
spare  coaxially  shielded  lead.    Such  a  line  is  constructed  by  passing  a 
stainless  steel  capillary  inside  a  0.6-cm-diameter  stainless  steel  tube 
extending  from  above  the  top  cryostat  flange  to  the  vacuum  seal  flange 
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above  the  He    evaporator.    Nylon  spacers  center  the  inner  capillary  which 
extends  through  a'  small  glass-to-metal  seal  mounted  in  the  bottom  end  of 
the  tube.    The  capillary  terminates  a  few  inches  below  the  top  end  of  the 
tube.    A  length  of  copper  wire  soldered  to  the  capillary  is  bent  in  a 
tight  spiral  before  extending  through  another  seal  mounted  in  the  upper 
end  of  the  tube. 

Eighteen  leads  made  of  formvar- insula  ted  #3^  Midohm  wire  enter  the 
vacuum  chamber  through  the  vacuum  chamber  pumping  line.    The  leads  terminate 
at  the  top  in  two  vacuum-tight  nine-pin  male  connectors.    All  leads  are 
thermally  grounded  to  the  vacuum  seal  flange  and  the  He**  evaporator. 
Only  five  of  these  leads  are  used  and  are  thermally  grounded  to  the  He3 
evaporator:    four  leads  serve  two  heater  resistors  and  one  lead  is  connect- 
ed to  a  nominally  U70-ohm  Speer  carbon  resistor  used  as  an  uncalibrated 
thermometer  during  cool-down. 

The  reduction  of  vibration  is  crucial  to  the  success  of  the  experiment. 
The  boiling  of  liquid  nitrogen  in  the  outer  dewar  (about  9  1.)  is  pumped 
continuously  with  a  Kinney  model  KC-k6  vacuum  pump  so  the  nitrogen  is  com- 
pletely solidified.    The  pump  is  capable  of  hO  cfm,  although  our  plumbing 
arrangements  for  pumping  on  the  dewar  are  awkward  and  reduce  this  figure. 
Nevertheless,  the  nitrogen  can  be  pumped  below  the  triple  point  in  about 
an  hour. 

The  cryostat  is  supported  entirely  by  a  partly  inflated  automobile 
inner  tube.    Considerable  vibration  can  be  transmitted  by  the  vacuum 
connections,  however,  unless  special  precautions  are  taken.    Vacuum  con- 
nection to  the  He1*  evaporator  is  made  through  a  four- foot- long  surgical 
rubber  hose  with  a  3/Vinch-i. d.    A  weight  suspended  from  the  ceiling  is 
clamped  to  the  hose  near  its  midpoint.     The  pumping    line  for  the  vacuum 
chamber  needs  to  be  quite  large.    After  this  chamber  has  been  evacuated,  a 
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valve  at  the  top  of  the  cryostat  is  closed  and  a  quick-disconnect  coupling 
in  the  vacuum  line  broken.    The  connections  for  the  sample  chamber  capillary, 
the  vapor  pressure  volume,  and  the  He^  gas  return  line  are  made  via  two-foot- 
lengths  of  3/8-inch-i.d.  stainless  steel  flexible  tubing.    The  He^  diffusion 
pump  is  attached  near  the  top  of  the  cryostat  and  is  supported  by  the  inner 
tube.    A  four-foot-long  section  of  1-inch-i.d.  stainless  steel  flexible 
tubing  connects  the  pump  to  a  flange  bolted  to  the  wall.    Another  three- 
foot-long  flexible  tube  connects  this  flange  with  the  backing  pump.  A 
weight  suspended  by  a  wire  from  the  ceiling  is  clamped  to  the  tube  near  its 
midpoint.    As  this  pump  vibrates  badly,  the  return  line  communicates  with 
the  plumbing  on  the  cryostat  support  frame  via  a  length  of  3/8- inch-diam- 
eter stainless  steel  flexible  tubing  which  also  connects  to  a  flange 
bolted  to  the  wall. 

Temperature  Measurement  and  Control 

Temperatures  were  measured  and  regulated  by  means  of  a  germanium  re- 
sistance thermometer  used  with  the  three-terminal  ac  bridge  shown  in  Figure 
3.    The  "three- terminal"  designation  is  appropriate  in  that  neither  lead 
of  the  thermometer  is  grounded,  and  ground  is  used  as  a  guard  point.  Re- 
sistances to  ground  do  not  greatly  shift  the  null  point  but  reduce  the 
sensitivity  to  unbalance.     This  arrangement  is  analogous  to  the  three- 
terminal  capacitance  bridge  described  elsewhere.    Resistance  of  the  leads 
is  included,  so  recalibration  is  required  after  any  change  in  the  wiring 
which  could  change  the  lead  resistance.     The  thermometer  must  be  calibrated 
in  situ. 

We  have  used  a  General  Radio  model  1493  ratio  transformer  (seven 
decades  plus  slidewire)  driven  by  the.ltOO  Hz    oscillator  of  the  phase- sen- 
sitive detector,  isolated  by  a  Gertsch  ST-100B  shielded  transformer.  The 
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unbalance  is  amplified  by  a  Princeton  Applied  Research  model  Ct-k  preampli 
fier  in  the  "low-Z  double-ended"  mode.    A  five-stage  RC  filter  precedes 
final  amplification  with  a  Princeton  Applied  Research  model  120  phase-sen- 
sitive detector. 

The  standard  resistor  is  a  2  kilohm  wirewound  resistor  sealed  in 
a  copper  tube  with  epoxy  and  clamped  in  a  heavy  aluminum  block  for 
thermal  ballast.    Electrical  connections  are  made  via  foot-long  lengths 
of  shielded  leads  constructed  by  passing  thin  stainless  steel  capillary 
through  stainless  steel  tubes.     This  arrangement  was  designed  to  reduce 
the  heat  leak  to  the  resistor.    The  entire  assembly  is  encased  in 
Styrofoam  for  further  thermal  shielding.    No  drifting  of  the  resistance 
has  ever  been  observed. 

No  provision  is  made  for  balancing  the  out-of-phase  component.  A 
variable  capacitor  placed  across  the  thermometer  leads  has  been  found  to 
work  satisfactorily,  but  the  additional  sensitivity  this  affords  is  not 
necessary  for  the  present  experiment  where  the  limiting  factor  is  the 
vibration- induced  capacitance  fluctuations. 

The  potential  across  the  standard  resistor  and  thermometer  in  series 
is  0.1  volts,  so  at  1  K  the  power  dissipated  in  the  thermometer  is  about 
0.4  |i  W  and  decreases  as  the  temperature  decreases.    The  power  dissi- 
pated at  k  K  is  approximately  kO  u  W.    Doubling  the  bridge  voltage  does 
not  cause  heating. 

The  detector  output  drives  a  single-stage  common-base  amplifier  which 
controls  the  current  through  a  heater  resistor.     This  resistor  is  of  the 
deposited  film  variety  so  its  resistance  is  not  strongly  temperature- 
dependent.    As  the  thermometer  temperature  drops,  the  heater  current  is 
increased  so  the  sample  chamber  is  warmed,  and  vice  versa.    An  adjustment 
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of  the  sensitivity  of  the  heater  current,  I,  to  detector  output,  V,  is  pro- 
vided so  that 


where  x  is  the  fraction  of  the  detector  signal  sent  to  the  amplifier  and 
I  is  a  constant  background  current.    The  parameters  x  and  I  are  easily 


The  response  is  not  instantaneous  as  temperature  is  not  transmitted 
in  copper  with  an  infinite  velocity.    Therefore,  if  x  is  too  large,  the 
regulation  becomes  unstable  and  the  temperature  oscillates.    For  the  same 
reason,  it  is  difficult  to  use  the  regulator  if  the  heater  is  too  distant 
from  the  sensor.    On  the  other  hand,  temperature  gradients  may  exist  across 
the  sample  chamber  if  the  heater  is  adjacent  to  the  thermometer.  The 
regulator  output  may  be  switched  between  two  resistors.    One  is  located  at 
the  bottom  of  the  sample  chamber  near  the  thermometer  (refer  to  Figure  7). 
The  other  is  mounted  on  the  side  of  the  He    evaporator.    The  more  distant 
heater  was  used  in  all  measurements  in  the  critical  region,  although  no 
temperature  drifts  were  observed  under  any  conditions  when  the  heater 
nearer  the  sensor  was  used. 

We  were  able  to  measure  and  regulate  temperature  to  within  10  u  K,  al- 
though our  determination  of  the  absolute  temperature  was  not  this  accurate. 

5  it- 
Temperature  calibration  is  achieved  by  means  of  He    and  He  vapor- 
pressure  measurements.    Liquid  is  condensed  in  a  copper  container  with  a 
volume  of  approximately  1  cc  which  is  attached  to  the  He^  evaporator.  To 
reduce  the  quantity  of  gas  used,  small  (0.l6-cm-  and  0 . 32-cm-diameter) 
tubing  was  used  for  the  lower  sections  of  the  vapor  pressure  line.    At  very 
low  temperatures  where  thermomolecular  correction  is  important,  the 


xV  +  I    >  0 


otherwise 


adjusted. 
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14. 

level  of  the  liquid  He    bath  must  be  high  enough  to  cover  the  junction  to 
0.62  cm  tubing,  as  the  temperature  of  this  junction  must  be  known.     To  re- 
duce the  possibility  of  non-equilibrium  between  the  thermometer  resistor 
and  the  liquid,  no  electrical  heating  is  used  during  calibration  runs. 
Readings  are  taken  at  approximately  50  mK  intervals  except  at  the  lowest 

temperatures  where  the  interval  is  decreased. 
k 

The  He    calibration  extends  from  the  lambda  point  to  k.2  K.     The  I958 

k 

He    vapor  pressure  scale  is  used.    A  Q.K%  correction  to  our  mercury  mano- 
meter pressure  readings  must  be  applied  to  give  the  equivalent  0  C  pressures 

No  calibration  data  are  taken  below  the  lambda  point  of  He\ 
3 

The  He    calibration  extends  from  below  Q.k  K  to  just  above  the  lambda 

point  in  liquid  He  .    A  Texas  Instruments  model  ll+5  Precision  Pressure 

Gauge  with  a  type  1  Bourdon  capsule  designed  for  pressures  from  0  to  300 

Torr  is  used.    The  quartz  tube  in  this  capsule  shows  hysteresis  effects 

so  that  even  when  the  zero  is  readjusted  immediately  before  the  lowest 

pressure  measurements,  accuracy  may  not  be  better  than  a  few  microns. 

Pressures  are  corrected  for  the  thermomolecular  effect  using  the  tables  of 

Roberts  and  Sydoriak5^  for  a  0.3l8-cm-radius  tube  extending  from  the  helium 

bath  to  room  temperature.    This  correction  is  negligible  above  0.6  K.  The 
3 

T62  He    temPerature  scale  is  used  to  convert  vapor  pressures  to  temperature. 

The  temperature  scale  near  0.9  K  where  the  transition  occurs  is  not  a 
closed  subject.    Clearly,  small  errors  in  the  vapor  pressure  scale  could 
lead  to  large  errors  in  the  expansion  coefficient  where  the  slope  of  the 
calibration  curve  is  important.    Such  errors  are  not  likely  as  the  equa- 
tion      from  which  the         scale  is  generated  is  semiempirical  in  that  the 
temperature  dependence  is  predicted  theoretically  and  the  data  is  used  only 
to  find  coefficients.    Slowly  varying  errors  displace  the  critical  tempera- 
ture but  do  not  change  the  shape  of  the  expansion  coefficient  curve  signifi- 
cantly. 
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As  small  resistance  increments  must  be  measured  to  determine  pre- 
cisely the  shape  6f  the  capacitance-versus-temperature  curve,  a  procedure 
to  find  the  temperature  as  a  smooth  function  of  resistance  is  required. 
This  function  must  satisfy  two  criteria.    First,  we  cannot  permit  fluctua- 
tions in  the  temperature-versus-resistance  curve  such  as  might  cause  an 
irregular  slope  in  the  capacitance-versus-temperature  curve.    Second,  the 
function  must  agree  with  the  calibration  points. 

It  is  observed  that  the  resistance,  R,  is  a  function  of  the  temperature, 
T,  obtained  from  the  vapor  pressure  measurements  of  the  form 

with  only  slight  curvature.    We  thus  attempt  a  machine  fit  of  the  form 

i-1 


N 

T  =  Z  A 
1=1  1 


R-.518 


This  form  with  N  =  5  fits  the  temperature  with  an  RMS  deviation  less  than 
1  mK  above  0.6  K.     The  scatter  in  the  vapor  pressure  measurements  is 
greater  than  this  at  the  lowest  temperatures.     Since  there  are  but  five 
smoothly  varying  terms,  we  expect  that  the  calculated  temperature  will  be 
a  smooth  function  of  the  resistance.     In  practice,  the  program  must  initially 
calculate  the  resistance  from  the  ratio  transformer  reading,  X,  using  the 
relation 

R  =  2K  ohm  x  (1  -  X)/X. 
The  large  number  of  manipulations  involved  allows  round-off  error  to  de- 
teriorate the  quality  of  the  fit.     The  program,  diagrammed  in  Figure  k, 
therefore  calculates  corrections  to  the  parameters  A  .    The  result  of 
this  program  is  a  relation  for  the  temperature  as  a  function  of  the  ratio 
transformer  reading  which  satisfied  our  two  criteria. 
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Figure  h.    Flowchart  of  the  Thermometer-Resistor  Fitting  Program 
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'Measurement  of  Small  Length  Changes 

Although  the  thermal  expansion  anomaly  at  the  Curie  temperature  is 

quite  large,  the  expansion  coefficient  changes  rapidly  with  temperature, 

so  in  order  to  define  the  shape  of  the  curve  carefully,  one  must  measure 

the  length  change  over  a  small  temperature  interval.    These  resulting 

length  changes  are  quite  small.    For  example,  our  measurements  for 

CuKgCl^-  2H^0  indicate  that  when  the  reduced  temperature  is  lo"'  and  the 

reduced  temperature  interval  is  2x10    ,  the  relative  length  change 

A  la/lsC^10    .    Furthermore,  the  available  sample  may  be  small  so  the 

length  change  is  less  than  0.1  A*.     Recent  reviews  by  Collins  and  White^ 
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and  Swenson  indicate  the  sensitivities  of  various  techniques  employed 
to  measure  dilitation  at  reduced  temperatures.  In  approximate  order  of 
increasing  sensitivity  these  methods  are  the  following:"^ 

1.  X-ray  diffraction  determines  the  lattice  parameters  as  a  function 

of  temperature.    A  relative  length  resolution,  A  1  /l  ,  of  approximately 

s  s 

IjxlO*    is  obtained. 

2.  The  two-terminal  capacitance  technique  has  been  used  to  give  a 
length  resolution  of  lo'^cm. 

3.  The  optical  interferometer  resolves  a  length  change  of  about  lo'^cm. 
These  methods  are  not  sufficiently  sensitive  for  our  needs.  Three 

techniques  are  capable  of  sufficient  resolution:^9 

1.  The  differential  motion  of  two  grids, ^°  one  of  which  is  coupled 
to  the  sample  by  means  of  a  lever  arrangement  to  increase  the  sensitivity, 
is  sensed  by  two  photocells.    Sensitivity  on  the  order  of  0.L&  has  been 
achieved. 

2.  The  variable  differential  transformer^1  has  been  used  with  a 
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resolution  up  to  0.05  *  under  special  conditions.    This  is  perhaps  the 
most  sensitive  method  available  for  thermal  expansion  measurement.  Unfor- 
tunately, it  is  quite  difficult  to  operate  and  is  limited  to  liquid  helium 
temperatures  where  the  coil  resistance  is  low  and  the  geometry  does  not  change. 

3.    The  capacitance  technique  has  been  used  to  measure  (  I-5  )  and  the 
equation  of  state    '  5  as  well  as  the  thermal  expansion  coefficient  as 

pioneered  by  White.  k    The  technique  has  also  been  used  near  room  tempera- 
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tures.       A  resolution  of  less  than  0.1  A*  is  possible. 

Capacitance  Technique 

_  65 

Thompson      suggested  the  use  of  a  "capacitance  probe"  which  would  con- 
sist of  a  guard  ring  and  a  disk  parallel  to  another  surface.  Sample 

chambers  suitable  for  low- temperature  differential  or  absolute  thermal  ex- 
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pansion  measurements  are  described  by  White. 

It  is  necessary  to  find  the  relation  between  the  capacitance  and  the 
thermal  expansion  coefficient  of  the  sample  for  a  typical  geometry  such  as 
that  shown  schematically  in  Figure  5.     It  is  convenient  to  divide  the  pro- 
blem into  three  parts.    First,  we  find  the  relation  between  the  thermal  ex- 
pansion coefficient  of  the  sample  and  the  temperature  derivative  of  the  plate 
separation.    Next,  we  assume  parallel  plates  to  find  a  in  terms  of  the  re- 
ciprocal of  the  capacitance.    Finally,  we  consider  in  an  approximate  fashion 
the  effect  of  nonparallel  plates. 

To  undertake  the  first  part  of  the  problem,  assume  the  plate  and 
crystal  supports,  the  sample  chamber  frame,  and  the  capacitor  plates  are  all 
constructed  of  copper  having  a  thermal  expansion  coefficient  given  exactly 
by  °Cu'    ReferrinS  to  Figure  5,  it  is  clear  we  need  only  consider  the  expan- 
sion of  copper  over  a  length  1£  equal  to  the  sample  length,  ls>  plus  the 
plate'  separation,  g.    Therefore,  we  have 


•Figure  5.    Simplified  Drawing  of  the  Sample  Chamber 

A.  Guarded  (fixed)  plate 

B.  Guard  ring 

C.  Moveable  plate 

D.  Sample 

E.  Sample  support 
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Upon  taking  the  temperature  derivative,  we  have 


dlc  dls 


Cu 


or 
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dl 
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We  will  find  that  the  neat  cancellation  of  thermal  expansion  of  construc- 
tion materials  over  comparatively  great  lengths  is  probably  not  valid. 

Next,  the  relation  of  ^  to  the  capacitance  and  its  derivative  must 
be  determined.    Referring  again  to  Figure  5,  note  that  the  guarded  plate 
has  radius  r  and  the  guard  ring  gap  is  of  width  w.    The  other  plate  is 
assumed  parallel  and  at  a  distance  g.     The  capacitance,  C,  between  the 
plates  is  given  by  Maxwell      as  follows: 

€  2 

C  "      °  Hr      .        w(l  +  w/2r)/r 

8  1  +  0.22w/g  J 


22w/g    J*  [2] 


Define  1  by 
eq  3 


r         €    _  2 

C  =     o  nr 

1 

eq 


and  let 

a  =     2  CI  +   h  ) 
r  2r  ' 

and  b  =  0.22w. 

Now  Eq.  [2]  bee  omes 


or 


IS 

leq_1  =  g"1  [1  +  a(l  +  b/g)"1], 


2g  =  (a  +  l)leq  -  b  +  {[(a  +  l)leq  -  b]2  +  hbleq)U2.  [3] 


We  assume  the  guard  ring  and  the  plates  have  the  same  expansion  coeffi- 
cient.   Then  we  find 

dT(^  =    dT  [(W  +  r)/r  *  l3  "  0- 
It  therefore  follows  that 

dl  =  0- 

de 

The  derivative  ~  is  found  directly  from  Eq.  [3].    For  our  expansion 
chamber,  the  values  of  w  and  r  are  given  by 

w  =  0.051  mm 

and  r  =  0.70  cm, 


so  that  a  =  1+xlO    .    Before  writing  down  the  expression  for  ~.  we  make 

dT 


.3  de 
the  approximation 

a  +  l^x-<l 

and  g/lgq  =  1  +  a/(l  +  b/g)^l. 

Then  we  have 

♦  1/2  [(1  +  b/leq)/(l  -  b/leq)  -  1]  %  [k] 

with 

^eq  2  dC~* 

dT    "  21eqaCu  +  £o  nr    W  t  5] 


During  the  experiment  described  here, 

C  =  20.6  pf 

and  1    =0.58  cm. 

s 

Therefore,  we  have 

b/1     =  0.17 
eq 

and  1  „  ■  0.0065  cm. 

eq 

Substituted  into  Eqs.  [k]  and  [5],  this  gives 

1    dg      1    dleq      ,  -h 

I    dl=  1    dT*    +  6x10  aCu 


with 


isdTq  "       °bu  +  0.36  Pf  |£ 

Finally,  substituting  these  results  into  Eq.  [l],  we  have 


a  =  1.017  ar    -  T  ^ 
Cu      1  dT 
s 


dC  1 


=  1.017  Qfcu  -  [(0.034  otu  +  0.36  pf        )  +  6xl0'4  Q^J 

=  0.98  ol    -  0.36  pf  — -  . 
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Now  we  come  to  the  final  part  of  the  problem,  the  effect  of  non- 
parallel  plates  on  -|.    For  the  purpose  of  this  discussion  neglect  fring- 
ing effects  so  the  infinite  parallel  plate  capacitor  formula  can  be  used. 
Suppose  the  upper  plate  is  tilted  about  the  y-axis.    The  x-axis  is  assumed 
to  lie  in  the  plane  of  this  plate  with  the  origin  of  the  coordinate  system 
at  the  center  of  the  plate.     Then  consider  the  area  2(r2  -  x2)1/2dx  at  a 
distance  (g  +  x  sin  e)  from  the  lower  plate,  where  e  is  the  angle  between 
the  plates  and  g  is  the  separation  of  the  centers.     The  projection,  dA,  of 
this  area  on  a  plane  parallel  to  the  lower  plate  is 
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dA  =  2(r2  -.x2)1/2  cos  0  dx. 


The  capacitance  in  the  present  approximation  is  given  by 
C  =  €o  [2(r2  -  x2)1/2  cos  e  ]/(g  +  sin  9)  dx 


=  2  n  €o  cot  e  -  [(g/sin  e  )2  -  r2]  1/2  }, 


c  sm2fl  gr2  eo  c°3  6 

8~4n€cose    +  c 
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Hence,  we  have 


nr2  €. 


s  s 


where 


#  =  C  tan  0 
2  Hr 

o 


If  we  define  1     as  before,  and  let 
eq 


A  ■  2r  sin  e, 


...  1  A 

then  ♦  -  f  ?  • 

eq 

We  cannot  evaluate  the  angle  0,  but  since  it  is  not  expected  to  change 

for  the  small  length  changes  observed  below  k  K,  and  since  1      (or  C)  is 

eq 

approximately  constant  in  this  temperature  range,  we  may  consider  that  the 
2 

factor  (1  -  $  )  cos  0  in  Eq.  [7]  introduces  a  small  but  constant  uncertainty 
in  the  determination  of  the  expansion  coefficient.     It  does  not  affect  the 
shape  of  the  expansion  curve,  however. 


To  estimate  an  order-of-magnitude  value  for  the  correction,  assume 
that  the  plates  are  set  at  room  temperature  with  an  error  A/g  less  than 
lA,  and  that  $  does  not  increase  greatly  as  the  temperature  is  reduced. 
Then  the  constant  error  introduced  into  the  determination  of  the  linear 
expansion  coefficient  is  less  than  17.  when  C  =  20.6  pf . 

It  is  possible  to  measure  the  angle  0  if  the  guarded  plate  is  re- 
placed by  two  concentric  plates.    At  room  temperature  with  the  plates  par- 
allel, the  area  ratio  is  determined  from  the  ratio  of  the  capacitance  of 
the  inner  plate  with  the  outer  one  grounded  to  the  capacitance  of  the  outer 
plate  with  the  inner  one  grounded.    The  guard  ring  is  always  grounded.  Then 
at  low  temperature,  the  angle  can  be  calculated  from  the  ratio  of  the  capa- 
citances.   Such  an  arrangement  was  used  in  the  present  sample  chamber,  but 
imaginary  angles  were  indicated  ocassionally.     It  is  probable  that  the 
method  used  to  construct  the  plates  allowed  the  inner  plate  to  drop  slightly 
with  respect  to  the  other  upon  cooling. 

Two  methods  are  generally  used  to  measure  small  capacitances.  The 

capacitor  may  be  made  a  part  of  a  resonant  circuit  so  length  changes  are 
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measured  as  frequency  changes.        We  have  used  an  alternate  method  which 
involves  direct  measurement  of  the  three -terminal  capacitance  by  means  of  a 
bridge  circuit.    One  uses  a  three -terminal  technique  so  that  stray  capaci- 
tances between  the  plates  and  the  leads  to  ground  shunt  the  detector  but  do 
not  greatly  change  the  null  point.    The  variable  component  of  the  bridge  may 
be  the  transformer  and/or  the  capacitance.    We  have  used  the  General  Radio 
Type  l6l5-A  capacitance  bridge  complemented  by  a  l^OU-B  lOOpf  standard 
capacitor  and  driven  by  a  General  Radio  1311-A  audio  oscillator  operated  at 
5  KHz  and  10  volts.     The  General  Radio  1232-A  tuned  amplifier  and  null  de- 
tector is  used  in  the  untuned  mode  as  a  preamplifier  for  a  Princeton 
Applied  Research  HR-8  phase  sensitive  detector.    The  bridge  employs  a  ratio 


transformer  tapped  at  intervals  of  one-tenth,  and  eight  standard  capaci- 
tors, six  of  which  are  used  at  once.    A  simplified  diagram  of  the  bridge  is 
shown  in  Figure  6. 

One  deficiency  of  this  system  is  its  dependence  on  ambient  temperature 
The  external  standard  has  a  temperature  coefficient  of  2x.l0~^  pf/C,  while 
the  stability  of  the  bridge  is  limited  by  the  temperature  coefficient  of 
the  10  pf  standard  which  is  less  than  5xl0_^  pf/C.    The  external  standard 
is  therefore  housed  in  a  Styrofoam-packed  container.    The  bridge  itself  can 
not  be  well  isolated  since  controls  must  be  accessible,  and  levers  or  shaft 
constitute  a  heat  leak  which  is  difficult  to  reduce.    The  bridge  is  there- 
fore enclosed  in  a  Styrofoam  box  with  extensions  on  the  balance  levers  ex- 
tending through  slots  cut  in  a  heavy-gauge  aluminum  front  panel.    The  tem- 
perature of  another  aluminum  plate  inside  the  box  and  beneath  the  bridge 
is  regulated  by  the  amplified  out-of-balance  signal  from  a  bridge  circuit 
with  a  thermistor  in  one  arm.    Although  the  temperature  of  the  room  was 
coarsely  regulated,  a  small  drift  of  the  measured  capacitance  was  sometimes 
observed  which  could  be  correlated  with  room  temperature  changes.  Such 
drifts  were  not  of  sufficient  magnitude  to  change  the  shape  of  the  measured 
expansion  coefficient,  however. 

The  expansion  chamber  used  in  the  measurements  reported  here  is  shown 
in  Figure  7.    The  sample  bears  directly  on  the  under  side  of  a  flat  disc 
which  serves  as  one  of  the  capacitor  plates.     This  plate  is  held  under  ten- 
sion against  the  sample  by  a  thin  brass  bellows  which  exerts  a  small  verti- 
cal force.    The  bellows  ends  below  the  sample  so  the  platform  upon  which 
the  sample  rests  is  larger  than  the  inner  diameter  of  the  bellows  (approx- 
imately 1.1  cm).    Since  the  platform  and  the  sample  is  not  limited  in  size 
by  the  bellows,  wide  samples  may  be  used.    Wide  samples  are  desired  since 
the  support  for  the  capacitor  plate  must  be  as  steady  as  possible.  Three 
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Figure  7.    The  Sample  Chamber 

A.  Flange  for  indium  gasket 

B.  Vacuum-tight  electrical  connector 

C.  Differential  screw  for  fine  adjustment  of 

plate  separation 

D.  Screw  for  coarse  adjustment  of  plate 

separation 

E.  Spring 

F.  Guarded  (stationary)  capacitor  plate 

G.  Guard  ring 

H.  Moveable  capacitor  plate 

I.  Sample 

J.  Nylon  post 

K.  Thin-walled  brass  bellows 

L.  Heater  resistor 

M.  Germanium  resistance  thermometer 
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posts  connect  the  bellows  to  the  plate.    The  spaces  between  the  posts  permit 
easy  access  to  the  platform  upon  which  the  sample  rests  so  that  it  is  easy 
to  determine  that  the  sample  is  properly  positioned.    The  posts  are  made  of 
nylon  since  the  lower  plate  must  be  electrically  insulated  from  the  sample 
chamber.    The  bellows  and  nylon  posts  serve  to  maintain  tension  on  the  upper 
plate  but  they  do  not  affect  the  plate  spacing. 

The  method  of  constructing  the  plate-guard  ring  assembly  was  suggested 
by  Dr.  C.  A.  Swenson  and  is  due  to  Dr.  G.  K.  White.    The  plate  is  machined 
as  a  plug  with  a  taper  of  about  1%.    Two  layers  of  1  mil  mylar  are  wrapped 
around  the  plate  before  it  is  tapped  into  a  mating  hole  in  the  guard  ring. 
Finally,  both  flat  sides  are  machined  smooth,  considerable  care  being  re- 
quired to  avoid  leaving  copper  in  the  gap.    One  can  check  that  the  plate 
will  remain  secure  at  low  temperatures  by  immersing  the  assembly  in  liquid 
nitrogen. 

The  height  and  angle  of  this  plate  may  be  adjusted  over  a  wide  range  to 
facilitate  installation  of  samples  with  slightly  varying  sizes  and  shapes. 
Provision  is  also  made  for  setting  the  plate  with  a  "negative  gap"  at  room 
temperature  so  that  a  sample  with  an  expansion  coefficient  much  larger  than 
that  of  copper  can  still  be  used  with  a  small  gap  at  k  K,  as  required  for 
high  sensitivity.    The  mechanism  supporting  the  plate  is  shown  in  Figure  6. 
Three  screws  support  a  brass  ring  from  the  sample  chamber  frame.  These 
screws  provide  a  coarse  adjustment  of  the  plate  position.    Three  differen- 
tial screws  support  the  plate-guard  ring  assembly  below  the  ring.  The 
motion  of  the  k-h8  section  in  the  brass  ring  is  partly  cancelled  by  the 
motion  of  the  2-56  section  in  the  guard  ring  so  the  net  displacement  is 
equivalent  to  a  screw  with  536  threads  per  inch,  but  our  arrangement  is  much 
stronger.    Three  springs  push  down  on  the  guard  ring  so  the  coarse  adjusting 
screws  are  held  firmly  against  the  expansion  chamber  frame.    Motion  due  to 
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play  in  the  threads  is  reduced  by  lightly  greasing  the  screws  so  that  at 
low  temperatures,  where  the  grease  solidifies,  the  threads  are  held 
tightly. 

The  plates  are  adjusted  by  trial-and-error  procedures.    With  the 
sample  in  place,  the  coarse  adjusting  screws  are  used  to  bring  the  upper 
plate  close  to  the  lower  one.    Then  the  differential  screws  are  rotated 
until  the  desired  capacitance  is  reached  and  the  plates  appear  to  be 
parallel.    A  strong  light  behind  the  plates  facilitates  this  adjustment. 
We  have  observed  that  the  net  contraction  of  the  CuKgCl^^HgO  sample  is 
similar  to  that  of  copper,  so  the  plate  setting  does  not  change  greatly 
between  room  and  liquid  helium  temperatures.     This  has  not  been  true  with 
other  materials  we  have  studied,  so  several  room  temperature  gaps  must  be 
tried  before  a  suitable    k  K  gap  is  attained. 

Eq.  [6]  indicates  that  the  measured  expansion  is  given  by  O.56  pf/C  AC 
for  small  capacitance  changes.     The  bridge  detects  a  AC  of  lo"^  pf  pro- 
vided one  figure  can  be  interpolated  on  the  HR-8  detector.     This  gives  a 
length  resolution,  Al  /l  ,  of  8.5x10  *°  when  the  capacitance  is  20.6  pf. 
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Vibration  of  the  plates  caused  fluctuations  of  two  or  three  in  the  inter- 
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polated  figure  so  the  resolution  was  reduced  to  about  3x10  .  This 
corresponds  to  a  length  change  of  approximately  0.1  A*  for  our  small 
(O.58  cm)  sample. 


CHAPTER  IV 


EXPERIMENTAL  RESULTS  AND  DISCUSSION 

Analysis  of  the  Data 

This  chapter  describes  how  the  capacitance-versus-resistance  data  are 

i  dl 
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analyzed  to  give  the  linear  expansion  coefficient,  a  =  —   ^jj-  ,  in  the  a 
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direction  of  a  single  crystal  of  CuKgCl^-S^O  from  O.k  K  to  k.2  K.  Parti- 
cular attention  is  given  to  the  data  in  the  critical  region  near  Tc  and  in 
the  high- temperature  region  above  2  Tc-    From  data  in  these  two  regions, 
information  is  obtained  about  the  stress  dependence  of  the  exchange  para- 
meters and  the  behavior  of  the  spin  correlation  functions  discussed  in 
Chapter  II.     In  addition,  the  stress  derivative  of  the  transition  tempera- 
ture is  determined. 

The  raw  data  consist    of  approximately  600  capacitance  values  versus 
resistance-bridge  ratio-transformer  settings.    Two  sets  of  data  cover  the 
range  from  0-35  K  to  k. 2  K.    Several  data  sets  were  taken  in  the  critical 
region  both  warming  and  cooling.    As  each  pass  defines  the  transition 
point  more  precisely,  the  next  pass  covers  a  smaller  temperature  range  with 
a  greater  point  density.    The  temperature  interval  on  the  final  set  is 
approximately  100  uK,  giving  a  capacitance  interval  of  about  10  ^  pf,  the 
limit  of  our  resolution. 

The  first  step  in  the  analysis  is  the  same  for  all  data  points.    The  " 
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temperature  for  each  point  is  machine  calculated  as  described  in  the  pre-' 
vious  chapter.    As  we  have  seen,  the  derivative  of  the  inverse  of  the 
capacitance  is  required  (Eq.  [6]  of  Chapter  III).    Hence  it  was  decided 
to  work  with  the  inverse  of  the  capacitance  from  the  outset.    The  computer 
program  is  designed  to  convert  the  raw  data  to  reciprocal  capacitance- 
versus- temperature  points.     This  information  is  punched  on  cards  and  printed 
for  further  analysis. 

Physical  shocks  of  the  expansion  chamber  cause  the  capacitance  to 
shift  discontinuous ly.    Such  shifts  sometimes  occur  when  valves  on  the  cry- 
ostat  support  frame  are  turned,  when  heavy  objects  are  dropped  on  the  labor- 
atory floor,  and  when  liquid  nitrogen  or  liquid  helium  is  transferred. 
In  addition,  the  standards  in  the  capacitance  bridge  drift  slowly.     It  is 
therefore  necessary  to  shift  the  inverse  capacitance  of  each  data  set  by  a 
constant  to  bring  the  curves  into  coincidence.    These  constants  correspond 
to  shifts  of  less  than  0.1%  in  the  capacitance  for  data  taken  over  a 
period  of  five  days.    Further  analysis  is  based  on  the  shifted  reciprocal 
capacitance-versus-temperature  data.    A  few  of  these  points  are  shown  in 
Figure  8. 

At  temperatures  greater  than  1.5  K,  the  thermal  expansion  coefficient 
is  expected  to  vary  as  T-2.    As  the  expansion  is  quite  small  in  this  tem- 
perature range,  the  thermal  expansion  of  the  chamber  cannot  be  neglected. 
Hence  the  measured  expansion  given  by  Eq.  [6]  of  Chapter  III  is  expected 
to  obey  a  relation  of  the  form 

a-  o.98aCu  =  -0.56  Pf  ^ 


=  aT"2  +  bT  +  cT5, 


vhere  the  T  and  T5  terms  are  included  for  the  electronic  and  lattice  con- 
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tributions,  respectively,  of  copper.    We  integrate  the  form  predicted  for 
the  expansion  coefficient  to  find  a  form  for  the  reciprocal  capacitance. 
The  coefficients  in  the  resulting  equation  are  found  by  a  least-squares  fit 
to  the  data.    When  the  fitted  equation  is  differentiated,  the  result  is 

I08(a  -  0.98  aCu)  =  2.16T  -  0.508T5  +  289T-2. 
67 

Carr  and  Swenson     have  found  that  in  this  temperature  range, 

108  aQu  =  1.7xlO-2T    +  1()"5T5. 

This  correction  is  much  too  small  to  account  for  the  observed  T  and 
coefficients.     The  formula  we  have  used  is  wrong  for  two  reasons,  however. 
First,  the  derivation  of  Eq.  [1]  of  Chapter  III  required  the  cancellation 
of  copper  expansion  over  lengths  greater  than  5  cm.    This  copper  is  of 
commercial  grade,  not  OFHC,  and  different  parts  were  cut  from  stock  pur- 
chased at  different  times.    Second,  the  screws  supporting  the  upper  capaci- 
tor plate  are  made  of  brass  with  unknown  expansion  coefficient. 

Although  the  total  measured  thermal  expansion  coefficient  changes  from 

UOxlo"  K      to  10x10    K     between  1.5  K  and  k.l  K,  this  change  is  well  de- 

-2 

scribed  by  the  three-term  fit  which  has  the  form  expected  for  a  T 

magnetic  contribution  as  well  as  an  electronic  contribution  proportional  to 

T  and  a  lattice  contribution  proportional  to  T^.    While  it  is  not  possible 

to  predict  the  sample  chamber  behavior,  it  is  not  unreasonable  that  the  T 
3 

and  T    terms  are  sample  chamber  effects. 

_p 

A  fit  from  1.5  K  to  3  K  gives  nearly  the  same  coefficient  for  the  T 
term,  although  the  contribution  from  the  T  and  T^  terms  at  3  K  is  only 
about  20%  of  the  total  measured  expansion  coefficient  and  decreases  rapidly 
with  temperature.    At  2.5  K  the  contribution  is  less  than  the  experimental 
error.     Therefore,  we  find 
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a  «  T*2  1.5  K    <  T  <  2.5  K. 


Such  behavior  is  not  inconsistent  with  the  data  to  k.2  K  as  discussed 

above.    The  results  of  the  fit  are  summarized  in  Table  3* 

At  temperatures  below  2.5  K,  where  the  correction  for  sample  chamber 
effects  is  negligible,  Eq.  [6]  of  Chapter  III  gives 

dC"1 

a  »  -  0.56  pf 


dT  ' 

In  principle,  one  could  take  simple  differences  between  points  to  find 

a<^')   -  ■  0.36  Pf  <C'»:1.-IC"1  ■ 

The  scatter  in  the  data  makes  this  procedure  useless  for  closely  spaced 
points. 

Outside  the  critical  region  where  the  shape  of  the  capacitance  curve 
varies  slowly,  we  use  a  local  fit  procedure  which  fits  N  consecutive  data 
points  to  a  linear  form 

C'1=  Oi    +V  Ti<T<Ti+N-l' 


so  that     C^CT,      )  =  -  O.56  pf  x  A^ 


Values  of  N  used  range  from  four  to  twelve  depending  on  the  scatter  and 

the  density  of  points.     The  program  prints  the  a.  and  T      .  for  any  desired 

i  av,i  * 

sequence  of  N's.    Figure  9  shows  the  expansion  coefficient  versus  terapera- 
ture.    Above  2.5  K,  the  fitted  T     component  is  shown  as  a  broken  line. 

This  procedure  is  not  valid  quite  near  the  transition,  for  the  shape 
of  the  capacitance  curve  varies  rapidly.     If  there  is  no  scatter  in  the 
measurements,  a  straight  line  fit  over  a  range  of  temperatures  gives  a 
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TABLE  3.    Results  of  the  Least-Squares  Fits 


T  >  1.5  K 


-8  2 
a    =  (290  +  10)xlO  K/T 


-2  -1 

10     <  e  <  10 


logarithmic  law 

» 

a    =    [  -  (470  +  30)  log  €  +  190s]  xlo'V"1      T  <  T 

c 

a  =   [  -  (470  +  30)  log  €  -    50a]  xlo'V1     t  >  T 


a 


c 


power  law 

a  =  [86. 8b  €'1/(2'°  ±  +  xlo'V1      T  <  Tc 

a  =  [30.5b  €-1/(1-3  ±  °'*>  +  230b]  xlo-V1     T  >  T 

c 


'Function  of  the  coefficient  of  the  logarithmic  term.  This 
value  corresponds  to  the  coefficient  given. 

b 

Function  of  the  exponent.    This  value  corresponds  to  the 
exponent  given. 
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A  C'1 

derivative  similar  to  the  quotient  of  differences         —  for  a  somewhat 
smaller  range  A  T.    Suppose  the  expansion  coefficient  follows  a  logarith- 
mic law  so  that 

C*1  =  P  +  Q  e  +  R  log  e. 

Then  for  A  €  «  €  , 
av 


¥F    =    ^+^>|av+<1+    7^>  dVR1°S£)| 


av' 


where  average  values  are  evaluated  at  the  midpoint  of  the  reduced  tempera- 
ture interval.    The  second  term  on  the  right  dominates  in  the  critical 
region,  so  that  we  have 

A  C"1               A  €  .  dC"1  , 

— t  ~  =  (1  +            )  .„  near  T  . 

AT       v         6  dT  lav  c 

av  1 

-2  -5 

For  example,  when  €      =10      and  A  e  =  ^xlO    .  the  difference  between 
'  av  ' 

A  C m  dC  m 

Jr-j —  and  — —  is  kO  %.  Clearly,  another  procedure  must  be  found  for  use 
in  the  critical  region. 

Before  proceeding  with  a  machine  fit,  it  is  useful  to  estimate  T  . 
Suppose  the  expansion  coefficient  follows  a  logarithmic  law  so  that 

a  =  A  log  e  +  B. 
Then  we  have 

.Act/A  | 


|T(a  +  Aa)  —  T  |  -  lO^/A|T(a)  -  T  | 


A  plot  of  T(a  +  A  a)  versus  T(a)  for  some  convenient  A  a  over  the  tempera- 
ture range  for  which  the  logarithmic  dependence  holds  is  called  a  "Mangels- 

68 

dorf    retarded  function  plot."       The  resulting  straight  line  intersects 
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the  line  T(a  +  A  a)  =  T(a)  at  the  Curie  temperature.    This  plot  determines 
the  transition  temperature  and  checks  the  logarithmic  dependence.  Unfor- 
tunately, this  technique  is  only  useful  in  the  part  of  the  critical  region 
where  the  local  fitting  technique  can  be  used  to  find  the  expansion  co- 
efficient, as  smooth  values  for  a  are  required.    The  plot  shown  in  Figure 
-2  -1 

10  for  10     <  6  <  10      shows  the  failure  of  the  logarithmic  function  for 
T  >  T  ,  and  predicts  the  transition  temperature  with  an  expected  error  of 
about  +  1.0  mK. 

A  simple  perturbation  procedure  gives  the  critical  behavior  of  the 
expansion  coefficient.    A  form  which  exhibits  approximately  the  observed 
temperature  dependence  is  selected  so  the  deviations  will  be  small.  The 
integral  of  this  form  is  the  equation  to  which  the  data  are  fitted.  A 
considerable  simplification  to  the  fitting  procedure  is  effected  by  treat- 
ing Tfi  as  a  known  parameter  which  is  tried  at  various  values  in  a  small 
range  about  the  temperature  predicted  by  the  Mangelsdorf  plot.     The  output 

of  the  fit  program  consists  of  the  value  of  D  =  (C~*),.  -  (C-*)  ,  versus 

fit  obs 

reduced  temperature  for  each  point  and  the  coefficients  in  the  equation  for 
the  expansion  coefficient.    Then  since 

d(C_1) 

a=-o.36Pf  — 25-2*!   .  Gfu  ,  o.36  pf 

the  correction  to  the  fit  may  be  calculated  by  taking  the  temperature  de- 
rivative of  the  difference,  D,  which  will  be  smoothly  and  slowly  varying 
if  the  choice  of  function  to  be  fitted  is  good.     The  derivative  is  found 
by  taking  slopes  of  a  hand-smoothed  plot  of  D  versus  T. 

This  procedure  has  been  used  for  the  form  which  fitted  the  specific 
heat,2^  i.e., 


a  +  =  A  log  €  +  B+. 


6k 


o  m  o 

<X>  CO  CD 

o  6  6 


(>i)    (DV  +  Dli 
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- 

The  RMS  deviation  describing  the  capacitance  fit  is  a  function  of  the  weight 
ing  factor  applied  to  the  points.    Since  the  supercritical  region  is  large, 
points  more  distant  from  the  transition  must  be  weighted  more  heavily.  A 
good  fit  was  found  with  a  weight  proportional  to  the  reduced  temperature. 
However,  the  fit  is  not  a  strong  function  of  the  exact  nature  of  the  weight- 
ing factor.    The  RMS  deviation  is  not  a  strong  function  of  the  Curie  temper- 
ature chosen  as  points  quite  near  the  transition  are  given  little  weight. 
The  result  of  the  fit  is  given  in  Table  3  and  is  shown  as  the  broken 

straight  line  on  the  semilog  plot  of  Figure  11  for  T  <  T£  and  Figure  12 

-2 

for  T  >  Tc.    Points  derived  using  the  local  fit  are  shown  for  10     <  €  <  10 

The  continuous  line  shows  the  value  of  a  near  T    calculated  from  the  results 

c 

of  the  fit  as  described  above. 

A  fit  to  the  form 

-l/n+ 

W   "  +B± 

has  also  been  made.    The  n's  are  treated  as  known  parameters  in  the  fit,  and 

a  wide  range  of  values  have  been  tried.    As  the  power  law  gives  a  greater 

deviation  in  the  supercritical  region,  this  fit  is  a  stronger  function  of 

the  weighting  factor  than  the  logarithmic  fit.    Therefore,  only  points  in 
-2  -1 

the  range  10     <  €  <  10       are  considered,  each  with  an  equal  weight.  The 
RMS  deviation  is  not  a  strong  function  of  the  critical  temperature  since 
points  quite  near  the  transition  are  not  included  in  the  fit.    The  expan- 
sion coefficient  was  not  calculated  near  Tc  as  the  deviations,  D,  are  large 
and  the  results  should  be  less  accurate  than  values  calculated  from  the 
logarithmic  fit. 

The  results  of  the  fit  are  given  in  Table  3  and  are  shown  as  the 
broken  curve  on  the  semilog  plot  of  Figure  13  for  T  <  T    and  Figure  Ik  for 
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T  >  T^.    Points  derived  using  the  local  fit  are  shown  for  10      <  €  <  10 
The  smooth  line  is  the  value  of  a  calculated  from  the  logarithmic  fit. 

Comparison  with  the  Specific  Heat 

A  comparison  of  the  results  of  the  high- temperature  fit  given  in 

Table  3  with  the  high- temperature  specific  heat  results  of  Miedema,  et 
20 

al. ,      leads  to  two  conclusions. 

1.  For  temperatures  above  2  Tc>  the"  nearest  and  next-nearest  neigh- 
bor spin  correlation  functions  have  proportional  temperature  derivatives 
and/or  the  exchange  constants  have  the  same  stress  dependence.  The 
analysis  of  the  critical  region  results  suggests  the  latter  is  false  so 
that  the  spin  correlation  functions  must  have  proportional  temperature 
derivatives . 

19 

2.  The  Callen  and  Callen  cluster  expansion  results      indicate  the 
correlation  functions  have  equal  temperature  derivatives  so  the  shift  of 
the  Curie  point  with  stress  can  be  calculated  using  Eq.[5]  of  Chapter  II. 
We  find 


 c 


(l.kO  ±  0.09)xl0_11K  cm2/dyne, 


where  the  subscript  1  refers  to  the  a  direction. 

Perhaps  the  most  striking  feature  of  our  data  is  the  severe  rounding 

of  the  lambda  peak.    The  supercritical  region  extends  to  a  reduced  tempera 

-2 

ture,  es,  of  almost  10    .     The  discussion  of  Chapter  I  led  to  a  relation 

between  €g  and  the  number  of  magnetic  ions,  N,  in  a  volume  defined  by  the 

spin- fluctuation  mean  free  path  as  follows: 

N^€  "5-^l06. 
s   
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This  corresponds  to  a  mean  free  path  of  about  500  A*  extending  over  about 

100  unit  cells.    The  problem  is  to  predict  the  asymptotic  behavior  of  the 

expansion  coefficient  of  a  perfect  sample  from  the  available  data  which  are 

-2  -1 

limited  to  the  reduced  temperature  range  10      -  10      where  the  spin  coher- 
ence length  is  not  on  the  order  of  the  temperature-independent  mean  free  path. 
A  log-log  plot  of  specific  heat  versus  reduced  temperature  (Figure  5 

of  Ref.  20)  indicates  that  the  behavior  determined  in  the  reduced  tempera- 
-2  -1 

ture  interval  10      -  10      differs  from  that  determined  in  the  range  10  - 
-2 

10      if  one  assumes  a  power  law.     It  is  not  correct  to  extrapolate  the 

asymptotic  behavior  from  observations  in  the  more  distant  decade.  This 

indicates  that  the  power  law  is  not  appropriate  over  a  wide  temperature 

range.    On  the  other  hand,  a  semilog  plot  of  specific  heat  versus  reduced 

temperature  (Figure  k  of  Ref.  20)  indicates  the  asymptotic  behavior  may  be 

-?  -1 

extrapolated  from  the  reduced  temperature  interval  10      -  10      if  a  log- 
arithmic form  is  assumed.    The  thermodynamics  of  Chapter  II  indicates  the 
critical  behavior  of  the  thermal  expansion  coefficient  should  be  linearly 
related  to  the  specific  heat,  so  these  conclusions  are  expected  to  apply 
to  the  expansion  coefficient  also. 

The  Mangelsdorf  plot  (Figure  10)  and  comparison  of  the  fitted  expan- 
sion to  the  measured  values  (Figures  11  and  13)  suggest  the  behavior  in  the 
critical  region  below  the  transition  is  given  satisfactorily  by  the  logar- 
ithmic law.    Above  the  transition,  the  situation  is  less  clear,  but  the  ex- 
pansion coefficient  appears  to  fit  a  power  law  better  (Figure  12  and  1^). 
The  previous  discussion  warns  against  the  extrapolation  of  the  power  law 
determined  in  the  reduced  temperature  interval  10      -  10    .    Hence,  we  use 
the  logarithmic  form  for  the  purpose  of  making  comparisons  with  the  specific 
heat  in  the  asymptotic  limit. 


A  comparison  of  the  results  of  the  logarithmic  fit  given  in  Table  3 
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with  the  critical-region  specific  heat  measurements  of  Miedema,  et  al. , 
leads  to  two  conclusions: 

1.  In  the  critical  region,  the  specific  heat  is  not  proportional 
to  the  linear  expansion  coefficient.    Therefore,  the  nearest  and  second- 
nearest  neighbor  spin  correlation  functions  do  not  have  the  same  tempera- 
ture dependence  throughout  this  region,  and  the  exchange  constants  do  not 
have  the  same  stress  dependence. 

2.  The  shift  of  the  Curie  point  with  stress  can  be  calculated  using 
the  ratio  between  the  asymptotic  terms  in  the  expansion  coefficient  and 
the  specific  heat.    We  find 


dT 

c 

do, 


-  (1.53  ±  0.05)xl0"11K  cm2/dyne, 


'J/1 


where  the  subscript  refers  to  the  a  direction.    This  result  agrees  within 
experimental  error  with  the  value  found  above  from  different  considerations, 

Unfortunately,  there  is  not  sufficient  information  to  evaluate  the 
strain  dependence  of  the  transition  temperature.    The  moduli  of  compliance 
are  not  known.    Although  the  stress  dependence  of  the  Curie  point  is  the 
same  for  the  b  direction  as  the  a  direction,  we  cannot  determine  this 
quantity  for  the  c  direction  without  performing  another  series  of  measure- 
ments.   To  estimate  a  value  for  the  elastic  moduli,  assume  isotropic  be- 


havior so  that 


'ii 


hi 

So, 


~  K/3, 
/0j/i*T 


where  we  use  the  compressibility  of  alum  as  discussed  in  Chapter  II. 


Then  Eq.  [10]  of  Chapter  II  gives 
\ 

7, 

°J/1,T 

dT 

where  the  subscript  1  indicates  the  a  axis.     If  we  further  assume  all 

~  oe 

are  equal  for  i  =  1,...,3,  we  find  as  an  order-of -magnitude  value, 


d  log  T  . 
6    c  7/3. 

d  log  V 

Considering  the  assumptions  made,  this  comes  quite  close  to  the  value  -10/3 
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Bloch      has  found  for  many  systems  coupled  by  superexchange. 

Summary  of  Conclusions 

Wood  and  Dalton  have  compared  several  critical  properties  of 
CuK2Cl^*2H20  with  the  results  of  exact  series  and  Pade  approximant  extrapo- 
lations above  the  transition  and  non- interacting  spin-wave  theory  at  low 
temperatures.    They  find  that  the  salt  is  a  close  approximation  to  a  bcc 
ferromagnet  with  first-  and  second-nearest  neighbor  Heisenberg  interactions 
with  spin  1/2. 

The  thermodynamics  of  the  lambda  transition  shows  that  the  thermal  ex- 
pansion coefficient  is  linearly  related  to  the  specific  heat  in  the  critical 
region.    This  is  consistent  with  our  observations.    Furthermore,  the  shift 
of  the  Curie  point  with  stress  can  be  calculated  using  the  ratio  between 
the  asymptotic  terms  in  the  expansion  coefficient  and  the  specific  heat. 
We  find 

=  (1-53  +  0.05)xl0"11K  cm2/dyne 
for  stress  in  the  a  direction. 
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The  procedure  of  Callen  and  Callen      treats  all  single-ion  crystal 
field  effects  and  two-ion  interactions.    The  expansion  coefficient  is 
given  entirely  by  the  temperature  derivative  of  three  types  of  spin  corre- 
lation functions.    When  this  procedure  is  applied  to  CuKgCl^^l^O,  the  only 
contributing  correlation  functions  are  <  S^'  S^>  and  <  S^'       >,  the  nearest 
and  second-nearest  neighbor  isotropic  spin  correlation  functions,  respec- 
tively.   The  specific  heat  is  described  by  the  same  two  functions. 

In  the  critical  region,  the  specific  heat  is  not  proportional  to  the 
observed  linear  expansion  coefficient.    Therefore,  the  two  spin  correlation 
functions  do  not  have  the  same  temperature  derivative  and  the  two  exchange 
constants  do  not  have  the  same  stress  dependence.     The  stress  dependence  of 
the  exchange  constants  is  expected  to  be  only  a  very  weak  function  of  tem- 
perature, so  the  latter  result  should  hold  for  all  temperatures  below  k  K. 

Our  measurements  suggest  the  expansion  coefficient  is  proportional  to 
the  specific  heat  above  1.5  K.     Therefore,  the  two  exchange  constants  have 
the  same  stress  dependence  and/or  the  two  spin  correlation  functions  have 
proportional  temperature  derivatives.    Measurements  in  the  critical  region 
suggest  the  former  is  not  true,  so  this  result  supports  the  cluster  approx- 
imation method  applied  to  the  problem  of  first-  and  second-nearest  neighbor 
Heisenberg  interactions  with  spin  1/2  which  predicts  the  two  spin  correla- 
tion functions  have  equal  temperature  derivatives  above  the  critical  region. 

Furthermore,  assuming  equality  of  the  temperature  derivatives  of  the 
two  spin  correlation  functions,  one  can  calculate  the  stress  dependence  of 
the  transition  temperature  from  the  ratio  of  the  expansion  coefficient  to  the 
specific  heat. 

The  result  is  given  by 
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=    -  (1.1+0  +  0.09)xl0      K  cm  /dyne 


for  stress  in  the  a  direction.    This  is  equal  within  experimental  error 
to  the  value  found  in  the  critical  region  from  entirely  different  consid- 
erations. 
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